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ABSTRACT 


We  have  studied  'stimul  ited  Raman  scattering  (SRS)  and  nonlinear 
index  effects  in  liquids  composed  of  spherically  symmetric  molecules.  ’ 
We  ccncludejthat  although  SRS  can  be  observed  without  self-focusing, 
nonlinear  index  effects,  presumably  because  of  molecular  redistribution, 
are  present  and  should  lead  to  beam  instabilities  at  somewhat  higher 
powers. 

We  also  have  studied  SRS  in  hydrogen  gas  pumped  by  a  giant 
pulse  ruby  laser.  We  determined  the  stability  conditions  for  various 
cell  configurations  and  lengths.  We  designed  and  tested  an  Stokes 
oscillator  which  reliably  produced  up  to  200  kW  of  Stokes  power  in  a 
nearly  diffraction-limited  fundamental  mode.  An  H2  Stokes  amplifier 
with  over  30  dB  small  signal  gain  was  built.  While  tests  indicated  that 
an  amplifier  output  brightness  of  over  10^  W-cm*^  sr~l  could  be  ob¬ 
tained  with  the  input  provided  by  the  H2  Stokes  oscillator,  only  10^ 
W-cm"2  sr-l  was  achieved  because  of  difficulties  in  coupling  the  oscil¬ 
lator  output  into  the  amplifier.  These  difficulties  resulted  from  the 
particular  configuration  used  rather  than  any  fundamental  limitation. 

Finally,  we  have  made  a  theoretical  study  of  the  third  order 
polarizability  applied  to  hydrogen,  and  the  resulting  coupled  Stokes - 
anti-Stokes  waves.  The  study  included  steady  state  and  transient  be¬ 
havior,  a.s  well  as  some  possible  explanations  for  the  anomalous  Stokes 
output  from  hydrogen  gas. 
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I. 


INTRODUCTION 


In  this  program  we  have  developed  a  prototype  Raman  oscillator- 
amplifier  and  demonstrated  high  brightness  and  brightness  gain  with 
moderate  pump  energies.  We  have  shown  that  this  prototype  has  the 
potential  to  achieve  brightnesses  in  excess  of  10*6  W-  om"‘  sr'*  with 
a  5  J  ruby  pump  pulse,  and  have  found  no  limits  to  further  extensions. 

We  feel  that  we  were  approaching  the  achievement  of  this  brightness 
potential  at  the  end  of  the  contract  period.  The  work,  accomplished  can 
be  conveniently  divided  into  two  major  areas:  (1)  the  search  for  strongly 
Raman-active  liquids  which  are  not  subject  to  nonlinear  index  effects, 
particularly  self -focusing;  and  (2)  the  use  of  gaseous  hydrogen  as  a 
Raman  medium. 

The  search  for  Raman-active  liquids  without  nonlinear  effects 
centered  on  those  composed  of  spherically  symmetric  molecules.  This 
work  was  reported  in  the  Semiannual  Technical  Summary  Report  and  is 
included  below  in  Section  II.  Extensions  of  the  theory  reported  in  the 
Semiannual  Technical  Summary  Report  on  the  nonlinear  index  of 
symmetric-molecule  liquids  to  general  polarizations  are  given  in  Ap¬ 
pendices  I  and  II.  This  work  led  to  the  conclusion  that  although  perhaps 
an  order  of  magnitude  improvement  was  to  be  found  in  nonlinear  index 
effects  by  employ'ng  liquids  with  spherically  symmetric  molecules,  these 
effects  still  placed  an  unsatisfactory  limitation  on  the  operation  of  a  high 
brightness  Raman  device.  The  origin  of  these  effects  in  spherically 
symmetric  liquids  is  presumably  the  molecular  redistribution  phenome¬ 
non  analyzed  In  Appendices  I  and  II.  1 

The  second  area  of  work,  covering  the  use  of  gaseous  hydrogen 
as  a  Raman  medium,  is  described  in  Section  III.  Nonlinear  index  effects 
in  Ho  gas  are  negligibly  small. 2  Preliminary  work  in  H;  gas  used 
parallel -beam  pumping  and  cells  with  nominally  zero  feedback  to  study 
Stokes  output  as  a  iunction  of  pressure  and  ruby  laser  pump  intensity. 
The  results  of  these  experiments  resemble  those  which  would  be  expected 
if  about  1%  feedback  existed,  and  show  typical  oscillator  behavior^;  how¬ 
ever,  there  is  no  known  mechanism  which  could  provide  this  much  feed¬ 
back.  There  are  other  possible  (but  unlikely)  origins  of  striking  anoma¬ 
lous  behavior  which  do  not  postulate  any  unknown  mechanisms  and  which 
are  consistent  with  observations  to  date;  one  of  these  is  discussed  in 
Section  IV.  However,  further  experiments  would  be  required  to  verify 
these  speculations.  Aside  from  physics  problems,  our  preliminary  ex¬ 
periments  gave  useful  information  on  the  stability  conditions  for  pumping 
a  hydrogen  amplifier  cell. 

Following  the  preliminary  H^  experiments,  we  built  and  tested 
a  high  brightness  device  consisting  of  a  H2  Stokes  oscillator  amplifier 
pumped  by  a  giant-  pulse  ruby  oscillator  amplifier.  A  brightness  gain 
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of  20  dB  over  the  Stokes  oscillator  was  obtained,  giving  a  brightness  of 
1014  W-cm-'i  sr”l.  Strongly  saturated  Stokes  amplification  was  not 
achieved  because  of  technical  difficulties  which  could  not  be  eliminated 
during  the  contract  period.  However,  unlike  other  high  brightness  de¬ 
vices,  the  H2  Stokes  oscillator -amplifier  has  no  obvious  fundamental 
limitations  on  size  and  power  output. 

Section  IV  outlines  the  theory  of  the  entire  third-order  polariz¬ 
ability  which  is  the  leading  contributor  to  nonlinear  effects  in  fluids. 

The  Raman  susceptibility  and  many  other  terms  emerge.  The  coupled 
Stokes -anti -Stokes  wave  relations  that  result  from  employing  this  polar¬ 
izability  and  a  plane  wave  pump  field  in  Maxwell's  equations  are  reviewed. 
The  magnitudes  of  the  nonresonant  nonlinear  polarizability  terms  cannot 
be  estimated  from  meager  and  uncertain  Kerr  data  on  H2.  Guessing  a 
value,  we  demonstrate  how  the  anomalous  threshold  behavior  can  be 
accounted  for  qualitatively.  These  non-Raman  terms  would  not  affect 
a  well-designed  high-brightness  H2  Raman  amplifier,  but  they  could 
affect  the  coherence  of  the  low-output  feeding  oscillator. 
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II. 


EXPERIMENTS  WITH  SYMMETRIC  MOLECULE  LIQUIDS 


We  have  conducted  several  experiments  to  study  the  SRS  effect 
in  liquids  composed  of  spherically  symmetric  molecules.  Careful  meas¬ 
urements  revealed  no  evidence  of  se’^-focusing  in  SiBr4<  However,  with 
a  multimode  pump,  we  have  observed  evidence  of  a  nonlinear  index  in 
the  form  of  frequency  spreading  in  the  Stokes  radiation.  Threshold  meas¬ 
urements  for  symmetric  and  nonsymmetric  liquids,  combined  with  cross 
section  measurements,  have  provided  a  basis  for  comparison  of  theo¬ 
retically  predicted  and  observed  gains. 


A.  Search  for  Self -Focusing  in  SiBr4 

Our  initial  observations  of  the  laser  and  Raman  radiation  exiting 
from  a  10  cm  path  length  of  SiBr4  indicated  some  filamentary  structure. 
However,  careful  control  experiments  showedthat  this  was  probably  a 
result  of  a  combination  of  weak  higher  order  transverse  modes  in  the 
GPL,  and  experimental  difficulties  encountered  in  imaging  through  the 
2  m  spectrograph.  In  order  to  obtain  more  conclusive  results,  we  set 
up  the  apparatus  to  analyze  the  polarization  of  laser  and  Raman  radia¬ 
tion,  using  circularly  and  linearly  polarized  pumps.  4  This  technique 
has  proved  to  be  a  very  sensitive  test  for  self-focusing  in  such  liquids 
as  nitrobenzene.  4  The  radiation  exiting  from  the  cell  is  analyzed  for 
polarization  and  wavelength,  and  a  pronounced  depolarization  of  the 
radiation  is  observed  to  accompany  the  onset  of  self -focusing  (SRS 
threshold)  with  a  circularly  polarizod  pump.  We  detected  no  depolariza 
tion  in  SiBr4  with  this  apparatus,  although  the  theory  of  self -focusing 
(regardless  of  the  underlying  mechanism)  indicates  that  there  should  be 
strong  depolarization  if  self-focusing  occurs.  This  negative  result  held 
ever,  far  enough  above  threshold  to  achieve  substantial  conversion  to 
Stokes  power  and  to  excite  at  least  five  orders  of  Stokes  radiation.  It 
should  be  noted  that  if  a  1/2  mm  diameter  pump  by  beam  was  used,  the 
Stokes  radiation  often  appeared  in  a  small  filament  about  40  to  50  pi  in 
diameter.  Presumably  this  results  because  of  spatial  narrowing  of  the 
Stokes  caused  by  the  variation  of  intensity  across  the  pump  beam,  and 
it  demonstrates  the  large  gains  achieved  in  this  liquid.  Because  the 
spatial  variation  of  the  pump  affects  the  gain  exponentially,  the  am 
plified  Stokes  beam  falls  off  mo.e  rapidly  as  a  function  of  the  radius 
than  the  pump  beam.  For  example,  if  the  pump  beam  is  gaussian,  the 
gain  has  the  form  G  =  exp  (qL),  where  a  =  o0  exp  (••r^/ro^).  The 
amplified  Stokes  wave,  normalized  to  the  on-axis  amplitude,  is  then 


S(r)  =  exp[aQL  (exp  (  r2/ro2)  •  1)]. 
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This  function  falls  to  1/e  at  a  radius  ij  such  that  f| 

rl2/ro2  =  ^in(1  “  l/aoL).  | 

a  i 

. 

Since  a  L  =*  10  for  our  pump  intensity  in  SiBr^>  r\/r  0.3;  a  puinp 
diamete°  2rQ  =*  0.  2  mm  gives  a  beam  size  of  SO  p?  The  pump  beam 
is  derived  from  a  1.  5  mm  aperture  and  a  40  cm  focal  length  lens,  and  ^ 

has  structure  of  this  size.  A  pump  variation  a  =  aQ  (1  -  gives 

a  similar  reduction  in  beam  size  for  aQL  «  10.  -* 


B.  Spectral  Observations 

Using  the  2  m  spectrograph,  we  measured  a  S^kes  shift  of 
225  cm"'  in  SiBr^  giving  a  first  Stokes  wavelength  of  705.  3  mp.  This 
would  be  a  symmetric,  or  ''breathing, "  mode  of  vibration  of  the  SiJ3r4 
molecule.  With  tetrabutyl  tin  (TBT),  which  should  also  show  spherical 
symmetry,  we  observed  a  shift  of  2920  cm-',  o-x responding  to  a  Stokes 
wavelength  of  871.  1  mp.  This  indicates  that  a  C-H  vibration  has  the 
lowest  threshold  in  TBT.  Because  the  higher  order  Stokes  shifts  in 
TBT  fall  outside  the  range  of  the  instrument,  only  the  first  order  Stokes 
radiation  was  observed  in  this  material.  With  SiBr^,  up  to  six  orders 
of  Stokes  radiation  were  observed. 

Our  giant  pulse  laser  topically  operates  ;n  several  longitudinal 
modes.  A  Fabry-Perot  interierogram  typically  shows  two  to  four  modes, 
with  spacings  of  about  0.6  cm"',  Other  observations  with  the  high  dis¬ 
persion  grating  of  the  Z  m  spectrograph  typically  indicate  two  modes 
spaced  by  0.  6  to  0.  7  cm"  1.  There  is  usually  a  time  modulation  of  the 
pulse,  indicating  other,  more  closely  spaced  modes,  probably  adjacent 
longitudinal  modes  of  the  laser  cavity. 

It  has  been  well  established  that  frequency  rendering  is  observed 
in  the  Raman  emission  of  liquids  which  self -focus  when  a  multimode 
pump  is  used.®’  °  This  results  from  a  modulation  of  the  nonlinear  index 
at  harmonics  of  the  difference  frequencies  of  the  pump.?  We  have  ob¬ 
served  3uch  broadening  effects  in  SiBr^  and  TBT,  with  the  amount  of 
broadening  increasing  with  increasing  pump  power  above  threshold,  and 
also  increasing  with  tfie  Stokes  order  for  SiBr^.  We  have  observed  up 
to  S.  5  mp  broadening  in  the  fifth  order  Stokes  radiation  from  SiBr4  and 
up  to  7.  4  mp  broadening  in  the  first  Stokes  radiation  from  TBT.  This 
spreading  of  about  130  cm"'  is  many  times  the  laser  mode  spacing  of 
about  0.  7  cm'  *,  indicating  that  there  is  modulation  of  a  nonlinear  index 
in  these  materials.  How  ver,  at  least  in  SiBr4>  this  nonlinear  index 
dees  not  i  iad  to  depolarization  of  a  circularly  polarized  pump,  mr  best 
test  for  self -focusing. 
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c. 


Threshold  Measurements 


Another  good  indication  of  the  presence  of  self-focusing  is  the 
theoretical  comparison  of  predicted  gains  using  the  incident  laser 
ntensity  with  the  minimum  gain  needed  experimentally  to  observe 
Stokes  radiation.  For  this  purpose,  we  have  made  threshold  measure¬ 
ments  on  several  liquids.  These  will  be  compared  with  measured 
spontaneous  Raman  cross  sections  in  Section  II- D.  For  the  threshold 
measurements,  the  GPL  beam  was  incident  on  a  1.5  mm  aperture. 
Following  this,  a  40  cm  focal  length  lens  weakly  focused  the  pump  beam 
in  the  10  cm  cell.  A  magnesium  oxide  screen  behind  the  cell  scattered 
the  output  radiation  to  calibrated  Stokes  and  laser  detectors.  The  Stokes 
detection  system  was  capable  of  detecting  about  1  W  of  Stokes  power 
from  the  cell,  except  for  SiBr^,  where  the  minimum  detectable  lower 
was  somewhat  higher  because  of  the  filters  used.  With  this  arrange¬ 
ment  the  entire  cell  length  was  within  the  depth  of  focus  of  the  40  cm 
lens,  so  that  the  resulting  beam  (about  0.  5  mm  in  diameter)  ."as  ap¬ 
proximately  a  plane  wave  throughout  the  cell.  The  thresholds  measured 
with  this  apparatus  are  given  in  Table  I.  The  error  limits  indicate  the 
scatter  in  the  data,  and  do  net  include  calibration  and  other  systematic 
errors.  However,  these  errors,  which  should  be  less  than  20%,  should 
not  affect  the  relative  thresholds.  A  linearly  polarized  pump  beam  was 
used  for  all  of  these  measurements. 


TABLE  I 

Measured  Thresholds 


Material 

Threshold  Power,  MW 

Silicon  tetrabromide 

0.29  ±0.02 

Tetrabutyl  tin 

0.  29  ±0.  03 

Nitrobenzene 

0.  050  ±0.  005 

Benzene 

0.  14  ±0.  01 

Carbon  disulfide 

0.  024  ±0.  002 

Carbon  tetrachloride 

0.  65  ±0.  03 

These  powers  should  be  multiplied  by  500  to  give  ‘.he  threshold 
intensity  in  MW/crn^.  The  above  results  for  nitrobenzene  and  benzene 
can  be  compared  with  Wang's  results*}  by  taking  into  account  the  dif¬ 
ferent  beam  sizes  and  extrapolating  his  results  to  a  10  cm  cell  length. 
We  must  also  take  into  account  the  fact  that  our  beam  is  nongaussian  in 
profile.  We  have  measured  a  spot  size  in  the  focal  plane  of  the  40  cm 
lens  which  is  four  times  the  diffraction  limit  for  the  1,  5  mm  aperture 
used.  Therefore,  if  we  take  a  "characteristic  transverse  radius  of 
curvature''^  of  the  laser  intensity  of  one  quarter  of  the  beam  radius, 
we  find  a  predicted  threshold  of  0.  045  MW  for  nitrobenzene  and  0,  16 
MW  for  benzene.  Thus,  our  results  for  these  two  materials  are  in  good 
agreement  with  Wang's  results  and  the  theory  of  self -focusing. 


D.  Cross  Section  Measurements 


In  this  section  we  describe  the  peak  normal  Raman  cross  section 
measurements  and  compare  the  resulting  theoretical  gains  with  the  above 
threshold  measurements. 

1.  Discussion  of  Experiments  and  Data 

The  normal  Raman  scattering  data  were  taken  with  a  cw 
argon  laser  operating  at  514.  5  mp  and  a  double  monochromator  arrange¬ 
ment  with  photodetection.  The  relative  cross  sections  of  several  liquids 
were  measured  and  the  absolute  peak  cross  sections  calculated  using  our 
previous  result  of  1.  3  cm~2  cross  section  per  unit  volume  per  unit  wave¬ 
length  for  nitrobenzene.  9  The  results  are  given  in  Table  II. 


TABLE  II 

Raman  Scattering  Data 


Liquid 

Raman  Line,  ^ 
cm'  * 

Peak  Cross 
Section, 

cm  "2 

Gain, 

cm"  *  /  MW  /  cm^ 

Nitrobenzene 

1345 

1.3 

1.4  x  10-3 

Benzene 

992 

2.  4 

2.  52 

Carbon  Tetrachloride 

440 

0.82 

0.82 

Tin  Tetrachloride 

370 

2.  9 

2.  54 

Silicon  Tetrabromide 

225 

5.8 

4.23 
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Unfortunately,  no  spontaneous  Raman  lines  were  observed  in  TBT  with 
this  apparatus. 

The  accuracy  of  these  measurements  is  not  well  known,  but  is 
judged  to  be  about  30%.  The  major  limiting  factor  is  that  the  double 
monochromator  used  had  an  instrumental  width  of  about  0.  07  mp,  mak¬ 
ing  it  difficult  to  determine  linewidths  accurately. 

2.  Comparison  of  Theoretical  Gain  with  Threshola 
Measurements  ~ 


The  largest  initial  intensity  available  for  amplification 
in  the  Raman  cell  results  from  the  quantum  mechanical  zero  point  vibra¬ 
tions  in  the  electromagnetic  field  and  is  of  the  order  of  10“^  W/cm  , 

For  our  beam,  this  represents  2  x  10"®  W  power.  Since  our  detection 
system  was  capable  of  detecting  about  1  W  of  output  power,  the  required 
gain  is  about  5  x  10?.  For  a  double  pass  through  a  10  cm  cell,  this  im¬ 
plies  a  gain  of  0.  88  cm"*.  Table  III  gives  theoretical  gains  obtained  by 
multiplying  the  above  gains  (cm^/MW/cm^)  by  the  threshold  laser  in¬ 
tensity  observed  experimentally.  F  is  the  factor  required  to  increase 
the  laser  powe*  to  give  a  theoretical  gain  equal  to  the  gain  required  ex¬ 
perimentally  for  the  observation  of  Stokes  emission. 


TABLE  III 
Theoretical  Cains 


Liquid 

Theoretical  Gain 

F 

Nitrobenzene 

0.  035  cm"1 

25 

Benzene 

0.  x  '6 

5 

Carbon  tetrachloride 

0.  266 

3.  3 

Silicon  tetrabromide 

0.  ol 5 

1.4 

With  the  many  cumulative  errors  involved  in  this  comparison, 
the  results  appear  adequate  to  explain  the  observed  SRS  from  SiBr^  on 
the  basis  of  purely  exponential  gain.  The  results  for  nitrobenzene  are 
consistent  with  our  earlier  measurements, ^  and  reaffirm  the  necessity 
of  invoking  the  self  .  ocusing  process  to  explain  the  observed  SRS.  The 
present  results  for  benzene  and  carbon  tetrachloride  are  not  conclusive. 
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Conclusions 


Our  measurements  and  theoretical  expectations  indicate  that  even 
symmetric  molecule  liquids  will  not  be  usable  at  desired  power  densities 
without  beam  instabilities.  The  theory  of  Appendix  II  shows  that  the 
threshold  power  for  the  self-focusing  instability  in  such  liquids  can  be 
accurately  inferred  from  dc  Kerr  effect  measurements.  The  lowest 
Kerr  constant  known  to  us  is  that  of  tetrabutyl  tin,10  for  which  a  thresh¬ 
old  power  for  self -focusing  of  several  megawatts  is  predicted. 
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III. 


EXPERIMENTS  WITH  RAMAN  EFFECT  IN  HYDROGEN  GAS 


Our  work  with  H^  gas  consisted  of  preliminary  experiments  to 
study  the  stability  properties  with  our  giant  pulse  ruby  laser  pump,  fol¬ 
lowed  by  the  construction  and  operation  of  a  hydrogen  Raman  oscillator- 
amplifier  designed  to  achieve  a  high  brightness  output. 


A.  Hyd  )gen  Stability  Experiments 

We  have  performed  a  number  of  experiments  with  H2  cells 
pumped  by  a  parallel  giant  pulse  ruby  laser  beam.  These  experiments 
typically  show  an  oscillator  type  instability  in  the  Stokes  output  as  a 
function  of  pump  intensity.-^  The  objectives  were  to  delineate  the  de¬ 
pendence  of  this  behavior  on  H2  pressure  and  the  physical  configuration 
of  the  apparatus.  In  particular,  we  sought  to  eliminate  sources  of  feed¬ 
back  which  could  lead,  to  the  oscillator  behavior. 


In  a  50  cm  long  H2  cell  with  uncoated,  nearly  perpendicular 
windows,  we  did  not  detect  Stokes  emission  (at  a  level  of  approximately 
1  W)  over  a  range  of  ~  14  to  70  atm  of  H2.  The  maximum  single-pass 
gain  calculated  was  exp  (4).  The  gain  G  =  exp  (p  IL)  is  calculated  using 
the  pressure-independent  gain  constant  p  =  1.  5  x  10“^  cm/MW,^  the 
measured  length  of  the  cell,  and  the  observed  laser  pump  power  con¬ 
verted  to  pump  intensity.  The  conversion  factor  depends  on  the  degree 
of  transverse  mode  control  used.  The  actual  powev  density  of  the  .aser 
beam  at  the  location  of  the  Raman  cell  was  checked  experimentally.  A 
0.  17  mm  diameter  aperture  was  carefully  scanned  across  the  laser  beam 
in  steps  of  0.20  mm,  and  the  laser  output  power  was  compared  with  the 
power  measured  on  the  other  side  of  the  aoerture.  The  transverse  mode 
structure  of  the  laser  was  monitored  for  each  data  point.  The  scan  wzs 
made  twice,  once  with  moderately  successful  efforts  at  transverse  mode 
control  and  once  with  no  effort  to  control  the  transverse  mode.  The 
scatter  in  the  data  was  quite  severe  for  the  n^n-mode-selected  case, 
with  noticeable  modulation  on  the  transmitted  beam.  The  mode-selected 
data  showed  some  scatter  caused  by  variations  in  the  effectiveness  of  the 
transverse  mode  control  and  a  much  smaller  residual  scatter.  Systematic 
variations  and  variations  which  are  n^t  understood  appeared  in  both  cases, 
but  were  much  less  apparent  in  the  mode -selected  case.  The  full  width 
at  half  maximum  for  the  two  cases  (from  the  smoothed  curves)  were 
1.  3  ±0.  2  mm  and  2.  9  ±0.  5  mm.  Both  curves  appeared  gaussian  in  shape, 
although  some  asymmetries  were  evident.  The  relevant  power  density 
values  were  12  ±1  MW/cm^  per  MW  of  power  output  and  3  ±1  MW /cm ^ 
per  MW  of  power  output.  Typical  power  output  values  were  5  MW  (mode- 
sel  ed)  and  10  MW  (non-mode -selected).  These  then  correspond  to  a 
power  density  of  60  MW/cm-  (mode -selected)  and  30  MW/cm^  (non¬ 
mode -selected). 
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In  a  90  cm  long  cell  with  uncoated,  nearly  perpendicular  windows, 
we  observed  Stokes  emission  at  a  calculated  gain  of  about  exp  (6).  The 
Stokes  power  increased  very  rapidly  with  pump  intensity,  a  character¬ 
istic  of  oscillator  threshold.  There  was  some  indication  of  a  lower 
threshold  at  the  higher  H2  pressures,  above  10  to  15  atm.  This  is  in 
contrast  to  the  strong  dependence  on  pressure  observed  by  other 

workers. 12,  1 3 

"We  attempted  to  reduce  feedback  by  tilting  the  cell  windows 
farther  from  alignment  with  the  pump  beam.  A  slight  increase  in  thresh¬ 
old  with  this  configuration  led  us  to  insert  a  dispersing  prism  between 
the  pump  laser  and  H2  cell  (with  tilted  windows)  to  block  feedback  of  the 
Stokes  radiation  through  the  pump  laser.  There  was  no  further  increase 
in  threshold. 

The  apparent  effect  of  tilting  the  windows  on  the  90  cm  cell  led 
to  the  hypothesis  of  "diffraction  feedback"  from  cell  windows  which  are 
not  far  from  alignment  with  the  pump  beam.  In  this  effect,  the  diffrac¬ 
tion  from  the  4%  reflection  of  the  spatially  narrowed  (see  Section  II-A) 
Stokes  beam  would  provide  the  source  of  feedback  leading  to  oscillator 
behavior.  Calculations  showed  that  under  typical  operating  conditions 
this  effect  gave  feedback  within  an  order  of  magnitude  of  that  required 
to  explain  the  oscillator  behavior.  In  order  to  check  this  hypothesis 
and  to  further  define  the  stability  criteria  for  our  H2  amplifier,  we  used 
a  150  cm  long  section  of  the  amplifier  cell  in  the  same  parallel  pump 
beam  configuration.  The  Bxewster's  angle  windows  on  this  cell  would 
reduce  the  diffraction  feedback  by  orders  of  magnitude.  However,  there 
was  little  effect  on  threshold,  with  detectable  Stokes  emission  appearing 
at  a  calculated  gain  of  about  exp  (7).  Thus,  while  diffraction  feedback 
appears  to  be  a  contributing  factor,  it  is  not  dominant  in  producing  the 
oscillator-like  instability  in  our  experiments. 

The  majority  of  our  preliminary  H2  experiments  were  performed 
with  the  150  cm  long  section  of  the  amplifier  cell.  With  a  more  sensitive 
Stokes  detection  system  we  were  able  to  observe  the  spontaneous  Raman 
scattering  and  check  its  proportionality  to  Hg  pressure.  A  plot  of  rela¬ 
tive  Stokes  output  power  as  a  function  of  ruby  laser  pump  power  is  shown 
.1  Fig.  1.  The  curves  show  a  linear  relationship  fitted  to  the  spontaneous 
scattering  data  and  the  calculated  exponential  behavior  of  the  stimulated 
scattering,  assuming  a  gain  constant  of  1.  5  x  10" 3  cm/MW  and  a  pump 
laser  intensity  of  10  MW/cm^  per  MW  of  pump  power.  The  actual  rela¬ 
tionship  between  pump  power  and  intensity  varied  from  shot  to  shot, 
from  about  5  to  12  MW/cm^  per  MW,  depending  on  the  purity  of  the 
transverse  mode,  and  resulted  in  the  scatter  of  the  experimental  points. 
However,  the  data  show  that  an  exponential  gain  constant  from  2  to  4 
times  larger  than  1.  5  x  10~3  cm/MW  is  required  to  give  the  observed 
dependence  on  laser  power,  if  amplifier  behavior  (i.  e.  ,  exponential 
buildup)  is  assumed.  The  more  likely  explanation  is  that  Fig.  1  demon¬ 
strates  the  extremely  rapid  (faster  than  exponential)  buildup  of  oscillator 
output  slightly  above  threshold. 
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Similar  measurements  at  21.4,  7.8,  and  5.  1  atm  pressure  with 
a  somewhat  more  reproducible  brightness  from  the  pump  resulted  in 
the  following  estimates  of  gain  constants,  assuming  exponential  buildup: 

p21  =  3,8x10  3  cm/MW 

p,  „  =  2.  4  x  l(f 3  cm/MW 
B c  ,  =  2.  6  x  10”3  cm/MW  . 

These  values  are  estimated  to  be  accurate  within  50%.  The  apparent 
higher  gain  at  21.4  atm  could  be  caused  by  the  higher  effective  pulsed 
gain  with  larger  Raman  linewidth.*  *  The  output  at  a  given  pump  level 
increases  with  increasing  pressure,  in  agreement  with  the  pressure 
dependence  of  the  spontaneous  scattering.  Thus  it  does  not  appear  that 
there  is  a  large  qualitative  change  in  the  Stokes  scattering  as  the  pres¬ 
sure  is  increased  above  about  10  atm.^ 

The  following  three  conclusions  can  be  drawn  from  these  pre¬ 
liminary  experiments:  (1)  the  parallel  beam -pumped  geometry  shows 
an  oscillator -like  instability;  (2)  there  is  a  range  of  stable  operation  in 
the  unsaturated  amplification  up  to  a  total  gain  of  approximately  exp  (6); 
and  (3)  this  behavior  is  not  strongly  dependent  on  H2  pressure. 


B.  Hydrogen  Oscillator-Amplifier  Experiments 

Following  our  stability  experiments  we  developed  a  hydrogen 
oscillator-amplifier  arrangement,  pumped  by  a  ruby  giant  pulse 
oscillator-amplifier.  The  objective  was  to  obtain  a  high  brightness 
output  from  the  amplifier  by  (1)  using  the  output  of  a  diffraction-limited, 
transverse  mode -controlled  Stokes  oscillator  as  the  amplifier  input; 
and  (2)  operating  the  amplifier  with  strongly  saturated  gain.  A  small, 
relatively  low  output  power  oscillator  can  provide  a  high-quality  output 
by  the  use  of  a  high-Q  cavity  with  some  means  of  mode  selection.  The 
oscillator  thus  provides  a  high  quality,  large  signal  input  to  the  ampli¬ 
fier.  The  amplifier  operates  in  such  a  way  that  the  effects  of  hot  spots 
in  the  ruby  pump  beam  are  minimized,  since  the  gain  increases  linearly 
with  pump  intensity  rather  than  exponentially  m  the  large  signal  case. 
With  this  arrangement,  we  have  achieved  20  dB  brightness  gain  in  the 
amplifier,  and  an  output  brightness  of  1  \  10*A  W  cm sr“*.  Because 
of  unanticipated  difficulties  encountered  in  coupling  the  H2  oscillator 
output  into  the  amplifier,  we  did  not  achieve  Hicngiy  saturated  ampli¬ 
fication  in  the  amplifier.  This  limited  the  hr.elo.e.ss  gain  through  both 
power  gain  and  far  field  divergence  angle  Our  work  has  had  the 
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following  results:  (1)  it  has  demonstrated  that  this  approach  to  obtain¬ 
ing  high  brightness  is  feasible;  (2)  it  has  shown  the  areas  where  further 
work  is  needed  to  optimize  this  approach;  and  (3)  it  has  indicated  how 
this  device  can  be  scaled  up  to  obtain  very  high  brightness. 

In  the  following  discussion  we  shall  first  describe  the  configura¬ 
tion  and  operating  parameters  of  the  apparatus,  and  then  discuss  in 
detail  the  characteristics  of  the  H~  oscillator  and  amplifier.  Finally, 
we  shall  discuss  the  results,  the  improvements  required  for  optimum 
operation,  and  the  possibilities  for  scaling  the  device. 

1.  Apparatus 

The  apparatus  is  shown  schematically  in  Fig.  2(a).  A 
rotating  prism,  Q-switched  ruby  laser  oscillator  (shown  in  Fig.  2(b) 
and  (c))produces  a  10  MW  of  power  in  a  40  nsec  pulse.  The  output  beam 
from  this  oscillator  is  enlarged  to  approximately  1  cm  diameter  by  an 
afocal  telsscope  and  is  then  amplified  in  a  12  in.  long  by  5/8  in.  diam¬ 
eter,  Brewster  angle  ruby  amplifier  crystal  (shown  in  Fig.  2(d))  pumped 
by  up  to  20  kj  of  energy  into  four  linear  flashlamps.  The  output  of  the 
ruby  amplifier  is  in  excess  of  200  MW  and  6  J  at  full  pump  energy. 
Approximately  2  to  3%  of  this  output  is  beam  split  off  and  used  to  pump 
the  H?  Stokes  oscillator  through  a  diaphragm  aperture  (not  shown  in 
Fig.  2)  and  a  beam-reducing  telescope  which  match  the  pump  beam  to 
the  TEMoo  mode  of  the  Ho  oscillator  cavity.  The  remainder  of  the  ruby 
amplifier  output  is  delayea  by  about  13  nsec  and  combined  with  the  H2 
oscillator  output  at  a  dichroic  mirror.  The  H2  oscillator  has  a  50  cm 
long  cavity  formed  by  a  10  m  radius  mirror  and  a  flat  mirror.  These 
mirrors  reflect  99+%  and  85%,  respectively,at  0.975  p,  and  transmit 
85%  and  70%,  respectively,  at  0.894  p.  The  H2  oscillator  output 
appears  in  a  beam  of  ~  1  mm  diameter,  with  a  typical  peak  power  of 
about  100  kW.  This  output  beam  is  expanded  in  a  10-power  afocal  tele¬ 
scope  to  match  the  size  of  the  ruby  amplifier  output  beam,  and  the  two 
beams  are  combined  at  a  dichroic  mirror.  The  combined  beam  is  en¬ 
larged  to  about  2.  5  cm  diameter  by  a  final  afocal  telescope  and  then 
enters  the  H2  Stokes  amplifier.  The  H2  amplifier  is  3  m  long,  with  a 
clear  aperture  of  about  3  cm.  The  cell  has  Brewster  angle  windows. 

The  two  H2  cells  are  connected  together  and  filled  with  pure  H2  gas  at 
pressures  up  to  50  atm. 

Detectors  monitor  the  power  output  of  the  ruby  amplifier  and  the 
power  input  and  output  of  the  H2  amplifier.  Calibration  of  these  de  - 
tectors  is  derived  from  an  ITT  FW  114A  Bi  Planar  photodiode  calibrated 
at  0.694  p  and  from  the  relative  response  of  the  S-l  photosurface  at 
0.  694  and  0.  975  p.  The  far  field  patterns  at  the  H?  amplifier  input  and 
output  are  photographed  with  a  high  quality  1.6  m  meal  length  telescope. 
The  brightness  gain  is  measured  by  attenuating  the  H2  amplifier  output 
far-field  monitor  beam  and  comparing  the  observed  far  field  patterns 
of  input  and  output. 
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Fig.  2(a)  Schematic  diagram  of  giant  pulse  ruby 
laser  oscillator -amplifier,  hydrogen 
Stokes  oscillator -amplifier  apparatus. 


Fig.  2(b)  Giant  pulse  ruby  laser  oscillator  head. 
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Fig.  2(c).  Giant  puls 
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Fig.  2(d).  Brewster  angle 
shown  with  refli 
not  shown. 
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2.  The  Hydrogen  Stokes  Oscillator 

The  H>  Stokes  oscillator  is  the  key  element  in  this  sys- 
te.n.  It  should  provide  a  diffraction-limited  beam  (TEMqq  mode)  wiih 
sufficient  power  to  saturate  the  H2  amplifier.  If  the  amplifier  has 
30  dB  small  signal  gain,  an  input  of  several  tens  of  kilowatts  is  required 
to  strongly  saturate  the  ruby  laser  pump  beam. 

For  transverse  mode  control  in  the  oscillator  we  used  a  rather 
high  Q  cavity  with  well-defined  modes, ^and  selected  the  fundamental 
(TEMqq)  mode  by  carefully  matching  the  ruby  pump  beam  to  that  mode. 
The  matching  is  obtained  by  an  aperture  and  a  oeam-reducing  telescope 
immediately  preceding  the  oscillator  cavity.  The  input  size  of  the  pump 
beam  is  varied  by  adjusting  the  power  of  the  telescope  and  the  size  of 
the  aperture.  The  output  size  of  the  pump  beam  is  controlled  primarily 
by  adjusting  the  focus  of  the  lelescope  and,  to  a  lesser  degree,  by  the 
aperture  size  and  telescope  power.  Pump  power  is  controlled  by  the 
aperture  size  and,  if  necessary,  by  attenuation.  With  proper  adjust¬ 
ment,  this  pumping  method  selectively  excited  the  TEMqq  mode  of  the 
cavity. 


If  pressure  broadening  is  dominant,  the  Stokes  gain  is  independ¬ 
ent  of  pressure,  while  the  linewidth  increases  linearly  with  pressure.^ 
If  the  pump  intensity  is  large  enough  at  a  g:ven  pressure,  the  Stokes 
buildup  rate  is  so  fast  that  some  Fourier  components  of  the  Stokes  signal 
lie  outside  th"  Stokes  linewidth,  leading  to  redu'ed  effective  gain.**  For 
normal  Stokes  Raman  scattering  (as  in  the  forward  direction,  where 
coupling  to  anti-Stokes  can  be  neglected),  th«.  peak  power  gsin  per  sec¬ 
ond  of  a  wave  packet  is  2ircAvV  where  Av(cm~l)  is  the  fui  ,.dth  of  the 
tr  ansition,  c  is  the  velocity  of  light,  and 
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here  2ircAvG  is  the  peak  gain  per  second  that  the  pump  would  produce 
if  the  linewidth  were  infinite  (see  eq.  (31),  Section  IV).  A  plot  of  V/G, 
i.  e. ,  the  ratio  of  effective  to  directly  calculated  gain,  as  a  function  of 
G  (Droportional  to  pump  power)  is  shown  in  Fig.  3.  The  arrow  shows 
the  point  at  which  the  calcula  ®d  power  gain  (cm"l)  equals  the  linewidth 
(cm*  i). 


At  pressures  of  approximately  20  atm  we  observed  that  the  pump 
power  required  to  reach  oscillation  threshold  was  much  larger  than  that 
calculated  using  the  known  cavity  losses  and  a  gain  constant  of  1.  5  x  10“^ 
cm/MW.  This  phenomenon  can  be  explained  by  the  effect  described 
above.  Thus,  establishment  of  an  oscillation  within  the  time  duration 
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Fig.  3.  A  plot  of  the  ratio  of  effective  to  directly  calculated 
gain  as  a  function  oj.  directly  calculated  gain  (pro¬ 
portional  to  pump  power).  At  the  etrrow,  calculated 
gain  equals  line  width. 
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of  the  pump  pulee  required  excessively  high  pump  powers  at  this  pres¬ 
sure.  V/e  therefore  studied  the  characteristics  of  the  oscillator  at 
several  pressures  up  to  70  atm.  We  found  that  at  -  50  atm  we  could 
obtain  reliable  Stokes  output  with  the  theoretically  predicted  threshold 
pump  intensity. 

The  calculated  passive  mode  diameter  for  our  oscillator  cavity 
is  1.  6  mm  (Ref.  14).  However,  we  observed  near-field  patterns  only 
about  •  mm  in  diameter.  The  far-field  patterns  also  corresponded  to 
a  1  mi  .  beam.  We  believe  that  this  effect  results  from  the  spatial  varia¬ 
tion  of  the  gain  in  the  cavity,  i.  e. ,  the  thape  of  the  pump  beam.  These 
effects  were  first  studied  by  Kogelnikl^  and  extended  to  spherical- 
mirror  cavities  by  Casperson  and  Yariv.  16  Calculations  of  the  mode 
size  to  be  expcc^d  with  our  operating  parameters  agree  very  well  with 
our  observations.  ^ 

A  further  effect  resulting  from  the  same  phenomenon  was  also 
noticed:  if  the  axis  of  the  pump  beam  and  the  axis  of  the  passive  cavity 
are  not  collinear,  the  pump  beam  can  geometrically  "pull"  the  axis  of 
the  mode  generated  in  the  cavity  away  from  the  passive  mode  axis.  This 
situation  is  shown  schematically  in  Fig.  4,  where  the  offset  is  caused 
by  an  angular  alignment  error  a  in  the  front  reflector,  leading  to  an 
offset  A  =  aR  between  the  passive  cavity  mode  and  the  pump  beam. 

Note  that  the  passive  mode  axis  must  always  be  a  radius  of  the  spherical 
mirror  and  must  be  normal  to  the  flat  mirror.  If  a  =  3  x  10“  ^  rad, 
then  A  =  0.3  mm  for  the  10  m  radius  mirror  we  used.  This  is  approxi¬ 
mately  our  alignment  accuracy,  and  gives  an  offset  of  one-third  of  the 
output  beam  diameter.  Another  cause  of  offsets  is  cylindrical  distor¬ 
tion  in  the  cavity  mirrors  caused  by  uneven  distribution  of  the  force  on 
the  mirrors  over  the  mirror  mounts;  we  have  observed  offsets  of  up  to 
3  mm  from  this  cause.  These  distortions  could  be  reduced  to  negligible 
values  by  careful  adjustment  of  the  mounts.  Regardless  of  the  cause  of 
offsets,  the  spatial  distribution  of  the  pump  causes  the  active  cavity 
mode  to  form  essentially  along  this  pump  axis  because  of  the  same  _ 

physical  process  studied  by  Kogelnik^  5  and  by  Casperson  and  Yariv.  ‘ 
Because  the  cavity  decay  time  is  longer  than  the  pump  pulse  length, 
much  of  the  Stokes  energy  built  up  in  the  cavity  while  the  pump  is  on 
decays  after  the  pump  shuts  off.  Furthermore,  since  the  pulling  effect 
disappears  with  the  pump,  this  energy  shifts  to  the  passive  cavity  mode 
which  is  offset. 

The  more  deleterious  of  these  two  effects  is  the  reduction  of  the 
mode  size,  since  in  this  case  a  higher  power  telescope  must  be  uned  to 
provide  the  amplifier  input.  This  in  turn  leads  to  more  sensitivity  to 
alignment  errors,  beam  offsets,  and  chromatic  aberrations  in  the  lenses. 
For  example,  a  0.  3  mm  offset  at  the  oscillator  output,  corresponding 
to  our  cavity  alignment  accuracy,  gives  a  7  mm  offset  at  the  amplifer. 
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and  additional  errors  exaggerate  this  result.  The  power  generated 
during  the  pump  pulse  is  pulled  according  to  the  above  discussion,  so 
that  the  offset  is  smaller.  E  'en  so,  vignetting  can  occur  fairly  easily. 

We  have  obtained  considerable  data  on  the  performance  of  the 
oscillator.  Much  of  this  work  used  the  giant  puise  ruby  laser  oscillator 
directly  as  a  pump,  instead  of  the  arrangement  shown  in  Fig.  2.  The 
only  significant  difference  in  the  results  is  that  the  oscillator  is  some¬ 
what  more  reliable  when  pumped  directly  by  the  GP  oscillator.  This 
apparently  results  primarily  from  the  accumulation  of  optical  damage 
to  the  beam  splitter  when  the  amplified  ruby  beam  is  used. 

The  evidence  that  a  nearly  diffraction  limited  (TEMqo)  mode  is 
generated  in  the  H2  oscillator  includes  beam  divergence  measurements 
and  near-field  patterns  on  Polaroid  film.  These  data  show  that  a  beam 
of  about  1  mm  diameter  is  formed  in  the  oscillator,  reduced  from  the 
1. 6  mm  passiv®  mode  size  by  the  effects  of  the  tightly  coupled  pump 
beam.  15-17  addition,  our  experience  has  shown  that  threshold  pump 
power,  the  time  buildup  and  decay  characteristics  of  the  Stokes  output, 
and  the  output  modulation  resulting  from  the  beating  of  adjacent  longi¬ 
tudinal  modes  of  the  cavity  are  all  strongly  correlated  with  the  character 
of  the  transverse  mode  generated.  A  mode  close  to  the  desired  TEMqo 
mode  is  associated  with  a  lower  threshold,  slow  buildup  and  decay  with 
a  single  time  constant,  and  a  stable  modulation.  A  poor  mode  usually 
occurs  with  a  relatively  high  threshold,  very  rapid  buildup,  and  faster 
decay,  usually  with  two  or  three  time  constants  apparently  involved. 

As  mentioned  earlier,  if  \,as  necessary  to  operate  with  Hz  pressures 
approaching  50  atm  in  order  to  achieve  reliable  generation  of  a  good 
mode.  Figure  5  snows  some  typical  Stokec  output  data.  Figure  5(a) 
shows  about  125  kW  peak  Stokes  power  from  a  cavity  formed  by  two  flat 
mirrors  with  the  same  reflectivities  as  the  mirrors  otherwise  used. 

(All  other  data  are  with  the  cavity  shown  in  Fig.  2.  )  The  decay  is  fairly 
rapid,  with  two  distinct  decay  rates  observed.  Figure  5(b)  and  (c)  show 
160  kW  and  100  kW  peak  powers,  respectively,  with  rapid  buildup  and 
decay,  and  multiple  decay  rates.  The  modes  associated  with  Fig.  5(a), 
(b),  and  (c)  all  had  poor  beam  divergence.  Figure  5(d)  shows  the  output 
power  associated  with  a  good  mode  which  had  small  beam  divergence. 
Peak  output  power  is  55  kW,  and  the  pump  power  is  only  one-third  to 
one-quarter  of  that  required  in  the  other  parts  of  Fig.  5.  Figure  6 
shows  the  output  associated  with  good  modes  on  a  more  expanded  time 
scale.  Peak  powers  are  about  50  kW.  The  slow  buildup  and  smooth 
decay  are  evident,  as  well  as  the  stable  modulation  resulting  from  mode 
beating.  The  decay  rate  agrees  with  that  calculated  from  the  cavity 
losses,  and  the  modulation  period  equals  the  calculated  cavity  roundtrip 
transit  time.  The  pressure  is  about  40  atm. 

Another  important  feature  in  the  H2  oscillator  is  a  time  delay 
between  the  peak  pump  power  and  the  peak  Stokes  output.  Figures  7  and 
8  show  the  ruby  pump  and  Stokes  output  as  measured  by  a  single  detector, 
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with  sensitivity  to  the  two  wavelengths  appropriately  adjusted.  The 
peak  Stokes  power  in  Fig.  8(c)  is  about  125  kW,  and  H2  pressure  is 
48.  5  atm.  It  is  apparent  that  the  time  delay  decreases  as  the  oscillator 
is  pumped  further  above  threshold,  reaching  a  minimum  value  of  about 
15  nsec.  Near  threshold,  the  oscillator  must  effectively  integrate  the 
entire  pump  pulse  in  order  to  build  the  spontaneous  Stokes  to  an  ob¬ 
servable  level,  so  that  the  delay  time  is  on  the  order  of  the  pump  pulse 
width.  Considerably  above  threshold,  the  gain  per  transit  of  the  cavity 
is  large,  so  that  a  fixed  number  of  transits  determines  the  time  delay. 
Because  the  oscillator  is  normally  operated  considerably  above  threshold, 
an  optical  time  delay  of  about  13  nsec  was  provided  for  the  H2  amplifier 
pump  beam,  so  that  the  peak  pump  and  signal  powers  would  coincide  in 
time. 


We  conclude  that  at  H2  pressures  between  40  and  50  atm,  this 
apparatus  can  produce  greate  *  than  100  kW  Stokes  power  in  a  near  dif¬ 
fraction  limited  1  mm  diameter  fundamental  mode  with  good  reliability. 
Threshold,  beam  divergence,  and  time  behavior  of  the  output  agree  with 
expected  characteristics,  although  the  mode  size  is  smaller  than  the 
calculated  passive  cavity  mode  size  because  of  the  spatial  variation  of 
the  pump  beam. 

3.  The  Hydrogen  Stokes  Amplifier 

The  calculated  unsaturated  gain  of  the  H2  amplifier  is 
exp  (6.  5),  or  28  dB  for  100  MW  of  input  pump  power,  assuming  a  uniform 
distribution  of  the  pump  power  over  the  aperture.  The  pump  beam 
actually  has  some  "hot  spots"  where  the  pump  intensity  is  higher  than 
the  average.  Since  the  unsaturated  gain  increases  exponentially  with 
pump  intensity,  these  hot  spots  would  lead  to  strong  nonuniformities  in 
the  amplified  beam  and  therefore  increase  the  beam  divergence  of  the 
output.  This  effect  can  be  avoided  to  some  degree  by  strongly  saturating 
the  amplifier  as  close  to  the  input  as  possible.  Operation  at  high  signal 
levels  also  provides  a  higher  conversion  of  pump  to  Stokes  radiation. 

We  operated  the  H2  amplifier  over  a  wide  range  of  pump  levels 
with  no  input  signal,  in  order  to  check  its  stability.  The  amplifier  was 
stable  throughout  the  expected  operating  range.  Even  with  calculated 
gains  of  over  60  dB  there  was  less  than  5%  conversion  to  Stokes  power. 
These  tests  were  made  with  the  ruby  amplifier  output  going  directly 
into  the  H?  amplifier. 

With  a  dichroic  mirror  in  place  we  first  tested  the  system  by 
eliminating  the  telescope  following  the  H2  oscillator  but  including  the 
2.  5-power  telescope  preceding  the  amplifier.  Thus  we  had  a  pump 
beam  in  the  amplifier  which  was  about  2.  5  cm  in  diameter,  while  the 
signal  beam  was  about  0.  25  cm  in  diameter.  We  achieved  brightness 
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gains  of  about  16  dB  and  output  brightness  in  the  range  0.4  to  1  x  10^ 

W  -  cm“2  8r-l.  With  the  10-power  telescope  added,  the  signal  beam 
diameter  was  increased  to  2.  5  cm.  We  observed  brightness  gains  of 
about  20  dB  and  output  brightness  of  about  1  x  10^  W  -  cm  ,jr*l. 

The  relatively  small  output  brightness  resulted  because  of  a  loss  of  input 
signal  strength  of  about  13  dB.  This  loss  will  be  discussed  below.  This 
meant  that  we  were  not  able  to  reach  high  conversion  efficiencies,  nor 
could  we  take  advantage  of  the  relative  smoothing  effect  of  strong  satura¬ 
tion  in  the  amplifier. 

The  total  available  operating  time  with  the  system  was  rather 
severely  limited  by  optical  damage  to  the  dichroic  mirror.  Four  units 
were  available,  all  of  which  had  "hard"  metal -metal  oxide  coatings. 

These  coatings  are  more  susceptible  to  giant  pulse  optical  damage  than 
some  other  types  of  coatings,  particularly  those  developed  specifically 
to  withstand  higher  pulsed  energies. 

We  have  considered  several  possible  causes  for  the  large  loss 
of  Stokes  signal.  First,  misalignment  of  one  or  both  telescopes  carry¬ 
ing  the  Stokes  signal  or  alignment  of  the  two  telescopes  to  different  axes 
could  have  caused  vignetting  of  the  beam.  It  should  be  noted  that  the 
alignment  telescope  used  to  align  the  system  had  to  be  reversed  in  its 
mount  for  the  alignment  of  both  telescopes;  this  involves  some  systematic 
error,  which  is  small  enough  to  be  unimportant  unless  compounded  by 
other  errors.  Second,  our  conception  of  the  H2  oscillator  mode  develop¬ 
ment  as  a  function  of  time  could  be  wrong,  leading  to  a  mismatch  of 
beam  sizes.  This  possibility  cannot  be  completely  eliminated  without 
making  time  resolved  studies  of  the  modes;  however,  we  have  consider¬ 
able  data  from  the  H2  oscillator  which  we  feel  are  consistent  with  our 
conception  of  its  operation.  This  includes  far-field  patterns,  near-field 
mode  photographs  and  burn  patterns,  time  resolved  power  output  meas¬ 
urements,  and  correlation  of  these  features  with  controlled  misalign¬ 
ment  of  the  oscillator  cavity.  Third,  Stokes  beam  offsets  can  cause 
vignetting  in  the  system.  The  two  primary  causes  of  these  offsets  are 
the  mode  shifting  associated  with  misalignment,  which  was  discussed 
above,  and  distortions  of  the  pump  beam  as  a  result  of  optical  damage 
at  the  beam  splitter.  We  have  observed  offsets  from  both  causes.  As 
discussed  earlier,  the  offset  associated  with  our  oscillator  alignment 
accuracy  is  about  one-third  of  a  mode  diameter.  Further  alignment 
errors  and  alignment  instabilities  increase  this  offset.  Major  distortions 
of  the  pump  beam  appear  to  cause  breakup  of  the  oscillator  mode,  often 
leading  to  two  lobe  modes  similar  to  a  TEMqi-  Finally,  we  have  con¬ 
sidered  chromatic  and  eff-axis  aberrations  in  the  telescopes,  leading  to 
spreading  of  the  beam,  as  a  possible  cause  of  loss  of  Stokes  signal.  The 
change  of  focal  length  of  lenses  with  wavelength  typically  changes  a 
telescope  which  is  afocal  in  the  visible  or  red  into  effect! vely  a  negative 
lens  in  the  near  infrared.  Lf  f  1  and  -f2  are  the  nominal  focal  lengths  of 
the  two  lenses  making  up  the  telescope,  the  resultant  effective  focal  length 
for  a  dispersion  An  =  -  6  is  given  by 
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where  n  is  the  visible  index  of  refraction.  For  a  typical  relative  dis¬ 
persion  of  6/(n  -  1)  =  0.  01  for  the  wavelength  change  of  interest,  our 
10-power  telescope  becomes  effectively  a  lens  with  -150  cm  focal  length. 
Off-axis  effects  would  be  most  important  at  the  negative  lens  of  the 
2.  5-power  telescope.  Since  this  lens  is  only  about  twice  the  size  of  the 
beam  and  is  located  about  2  m  from  the  10-power  telescope,  the  cumula¬ 
tive  effect  of  offsets  and  alignment  errors  can  be  quite  large. 

Figure  9  shows  some  near  field  photographs  of  the  Stokes  oscil¬ 
lator  output  beam.  The  top  row,  Fig.  9(a)  through  (d),  shows  the  beam 
with  the  10-power  telescope  removed  and  the  2.  5-power  telescope  in 
place,  and  indicates  the  progressive  breaking  up  of  the  beam  with  added 
misalignment  or  pump  beam  distortion.  The  beam  size  is  about  3mm. 

A  reduced- size  TEMqo  mode  should  decay  primarily  through  the  TEMqo> 
TEM()2>  and  TEM20  modes  of  the  passive  cavity  after  the  pump  is  re¬ 
moved;  this  feature  can  be  seen  in  Fig.  9(a)  through  (c).  The  pattern  in 
Fig.  9(d)  occurred  after  the  beamsplitter  had  been  damaged.  Since  no 
aperture  of  the  system  is  included  in  these  photographs,  offsets  from 
the  alignment  axis  cannot  be  deduced.  The  bottom  row,  Fig.  9(e)  through 
(h),  shows  the  beam  with  the  10-power  telescope  in  place  and  the  2.  5- 
power  telescope  removed.  These  photographs  are  overexposed  to  show 
the  aperture  of  a  beam-folding  prism  after  the  10-power  telescope,  and 
therefore  exaggerate  the  beam  size  and  edge  structure.  The  beam  size 
is  somewhat  larger  than  1  cm.  Ideally,  the  beam  should  be  centered  in 
the  prism  aperture.  Offsets,  distortion,  and  breakup  of  the  beam  can 
be  seen.  The  decaying  Stokes  radiation  should  overfill  the  prism  aper¬ 
ture,  and  presumably  provides  the  illumination  of  that  aperture  in  these 
photographs. 

With  either  the  10-power  telescope  or  the  2.  5-power  telescope 
alone,  the  Stokes  signal  is  not  significantly  attenuated.  With  both  tele¬ 
scopes  in  the  beam,  about  13  dB  attenuation  is  observed,  and  photo¬ 
graphs  of  the  beam  show  that  the  aperture  of  the  telescope  is  grossly 
overfilled  following  the  negative  lens.  Based  on  the  above  discussion, 
we  conclude  that  offsets  due  to  alignment  errors,  beam  spreading  due 
to  focusing  errors,  errors  in  aligning  a  system  of  two  telescopes 
separated  by  about  2  m  with  three  beam-folding  prisms  between  them, 
and  the  limited  apertures  of  some  components  combine  to  move  the  beam 
off  axis  and  spread  and  attenuate  it  to  the  degree  we  observe,  even 
though  no  error  is  particularly  damaging  in  itself. 
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4. 


Discussion 


The  H2  Stokes  oscillator -amplifier  apparat'  is  rather 
complicated,  in  the  sense  that  it  has  some  components  whose  alignment 
is  critical.  This  makes  the  system  dif-'icult  to  align  and  allows  the 
development  of  cumulative  errors. 

Several  modifications  can  be  made  to  our  system  in  order  to 
alleviate  the  problems  we  encountered.  These  include  the  following: 

a.  Increase  the  mode  size  of  the  H2  oscillator.  This  would 

considerably  simplify  coupling  of  the  oscillator  output 
into  the  amplifier.  For  a  mirror  spacing  d/2  which  is 
small  compared  with  the  radius  R  of  the  spherical  mir¬ 
ror,  the  mode  diameter 


increases  rather  slowly  with  d  and  P.  For  large  R, 
alignment  and  rtability  become  problems.  One  solution 
would  be  ‘o  us'  a  confocal  or  nearly  confocal  cavity. 

This  would  give  the  maximum  spot  size  for  a  given  radius 
mirror  and  very  greatly  reduce  alignment  and  stability 
problems.  In  addition,  the  output  could  be  taken  from 
the  spherical  mirror,  further  increasing  the  output  mode 
si.  e.  Another  approach  would  be  to  use  a  small  output 
mirror,  leading  to  a  small  F  **snel  number  and  relatively 
high  diffraction  loss.  Such  diffraction  loss  is  quite  effec¬ 
tive  in  selecting  the  fundamental  mode  and  the  diffraction 
from  the  fundamental  mode  can  provide  the  output  beam. 

b.  Increase  the  transmission  of  the  output  mirror  of  the  H2 
oscillator  to  match  the  cavity  decay  time  to  the  pump 
pulse  duration.  This  step  would  supply  all  of  the  energy 
stored  in  the  cavity  as  useful  power  to  be  amplified  by 
the  pump  pulse.  The  pump  power  to  the  oscillator  of 
course  would  have  to  be  increased  to  compensate  for  the 
added  loss. 

c.  Decrease  the  distance  betv  ;en  the  two  telescopes  carry¬ 
ing  the  Stokes  beam.  This  would  simplify  coupling  of 
the  H2  scillator  output  into  the  amplifier.  In  addition, 

the  rearrangement  would  be  made  so  that  the  entire  system 
could  be  aligned  without  reversing  the  alignment  telescope 
in  its  mount,  eliminating  that  source  of  error. 
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d.  Use  a  high  field  damage  resistant  dichroic  mirror.  In 
our  apparatus  this  mirror  is  a  flat  with  high  reflectivity 
in  the  red.  Such  mirrors  are  available  with  high  field 
damage  resistant  characteristics  as  a  standard  item. 

e.  Partially  diffuse  the  H2  oscillator  pump  beam.  This 
would  provide  a  more  uniform  and  easily  controllable 
pump.  Using  some  of  the  above  techniques  of  mode  con¬ 
trol,  the  pump  beam  need  not  be  so  constricted,  and  the 
effects  of  the  rapid  spatial  variation  of  the  gain  would  be 
circumvented. 

The  primary  advantage  of  this  design  principle  is  that  because 
a  gaseous  amplifying  medium  is  used  and  the  pump  is  not  introduced 
from  the  side,  the  system  can  be  scaled  up  to  larger  size.  The  only 
requirement  is  that  the  intensity  be  kept  below  the  level  where  the  am¬ 
plifier  becomes  unstable.  An  amplifier  chain  thus  could  be  constructed, 
separated  by  afocal  telescopes  to  progressively  increase  the  beam  size 
and  dichroic  beamsplitters  to  introduce  additional  pump  energy.  Thc 
amplifier  chain  could  include  backward  wave  amplifiers  in  the  final 
stages  to  yield  very  high  power  in  a  short  pulse. 
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IV.  THEORY  OF  THIRD  ORDER  POLARIZABILITY 
APPLIED  TO  HYDROGEN 


A.  Introduc  tion 

The  major  theoretical  problem  of  the  Hg  Raman  laser  which 
we  have  approached  is  that  of  understanding  a  rather  complicated 
anomalous  Raman  output  observed  in  parallel-beam-pumped  H2  gas 
cells  with  nonreflecting  tilted,  or  Brewster  angle,  windows.  The 
output  is  anomalous  in  the  sense  that  it  cannot  be  understood  on.  the 
basis  of  the  known  Raman  transition  strength  and  measured  cavity 
parameters.  The  anomalous  behavior  occurs  at  pressures  above 
~  10  atm,  is  insensitive  to  the  higher  pressures  (up  to  ~  100  atm), 
and  resembles  that  which  would  be  expected  if  reflecting  windows  (of 
around  1%  reflectivity)  were  present. 3,  11,  13  We  have  considered 
dozens  of  possible  sources  of  this  behavior,  from  the  effects  of  non¬ 
linear  absorption  in  impurities  in  the  H2  gas  to  stimulated  normal 
and  thermal  Rayieigh  wing  scattering.  Some  of  these  possibilities 
have  been  discussed  in  previous  reports.  Upon  further  investigation, 
we  have  been  able  to  discard  most  of  such  hypotheses.  One  rather 
complicated  hypothesis,  explained  in  Section  V-F,  must  be  investi¬ 
gated  further.  It  presupposes  the  existence  of  a  nonresonant  non¬ 
linearity,  for  which  there  is  little  evidence,  and  involves  dielectric 
guiding  effects  stimulated  by  the  pump  beam  at  Raman  shifted  fre¬ 
quencies.  Independent  measurements  of  electronic  nonlinearities  in 
H2  would  be  the  most  sensitive  test  that  could  be  devised  for  this 
hypothesis  at  present.  Such  an  effect  need  not  affect  H2  Raman  laser 
performance  if  beam  diameters,  pressures,  and  oscillator  param¬ 
eters  are  adjusted  properly. 

In  the  following  sections  we  outline  the  derivation  of  the  appro¬ 
priate  nonlinear  susceptibility  and  the  coupled  wave  theory  appropri¬ 
ate  to  weak  stimulated  Raman  scattering  from  a  plane  wave  purr.p. 

All  parameters  are  evaluated  for  H £  according  to  the  best  present 
knowledge.  The  role  of  the  nomesonant  nonlinearities  is  established, 
and  limitations  imposed  by  transient  effects  are  discussed.  We  begin 
this  rather  self-contained  view  of  the  theory  with  an  important  sim¬ 
plifying  theorem  on  the  form  of  the  interaction  Hamiltonian. 


B.  A  Simplification  in  the  Form  of  the  Field-Molecule 
Interaction 


From  the  most  fundamental  (nonrelativistic )  viewpoint,  the 
electromagnetic  field  couples  to  a  charge  in  a  moiecule  via  a  term 
in  the  Hamiltonian  of  the  familiar  form  (in  cgs  units) 
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where  e  and  m  are  the  charge  and  mass  of  the  particle,  p  is  its 
momentum  operator,  and  the  A(x,  t)  and  $  (x,  t)  are  the  vector  and 
scalar  potentials  of  the  electromagnetic  field  evaluated  at  the  position 
x  of  the  particle.  Calculation  of  the  polarization  to  third  order  in  A 
and  £  with  this  interaction  leads  to  a  rather  complicated  sum  of 
delicately  cancelling  terms  which  are  difficult  to  use  for  numerical 
calculations.  We  have  been  able  to  show  that  in  the  dipole  approxima¬ 
tion,  this  exact  method  of  calculating  the  nonlinear  polarization  is 
equivalent  to  a  much  simpler  one  in  which  the  above  Hamiltonian  is 
replaced  by  only  one  term  -  m-  E  (Rt)  where  m  is  the  dipole  moment 
operator  for  the  molecule  aivTgT&t)  is  the  electric  field  at  its  center 
position  R  .  This  equivalence  has  not  been  evident  to  some  workers. 
We  shall  employ  the  worthwhile  simplification  and  numerical  con¬ 
venience  of  this  latter  form. 


C.  A  Quantum  Mechanical  Expansion  of  the  Nonlinear  Polarization 
to  Third  Order  in  the  Perturbing  filectric  Fields 

We  seek  the  expected  value  M^(t)  of  the  i^1  component  of  the 
electric  dipole  moment  of  an  isolated  molecule  in  an  electric  field 
whose  components  are  E:(s).  If  the  interaction  Hamiltonian  is 
-  m^(t)E^(t)  in  the  unperturbed  Heisenberg  representation  (repeated 
indices  are  assumed  to  be  summed),  a  well-known  quantum  mechanical 
perturbation  procedure  gives 

t 

M  (t)  =  i  J  ds  <  m.(t),  m.(s)  >  E  (s) 

J  J 

-  t  8  U 

+  i  /  ds  J  du  j  dv  <  [  [  [m.(t),  m.(s)]  ,  mk(u)]  ,  m^  (v)]  > 
•  E^(s)Ek(u)Ei  (v)  f  .  .  .  ,  (1) 

in  which  energies  are  measured  in  frequency  units  (h  =  1 )  and  the 
brackets  <>  indicate  the  expected  value  for  an  unperturbed  (Ej  =  0) 
molecule.  The  lower  time  limits  on  the  integrals  are  some  time  in  the 
distant  past  before  which  the  fields  were  absent  and  the  molecule  had 
the  probability  P&  of  being  in  an  energy  eigenstate  |a>. 
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We  shall  assume  that  the  electron  field  is  a  superposition  of 
oscil’ating  components  ZpE^p  exp  iwpt,  of  (plus  and  minus)  frequencies 
wp,  and  complex  amplitudes  E^p,  where  t*ie  amplitude  of  a  minus  fre¬ 
quency  term  is  the  complex  conjugate  of  the  corresponding  positive 
frequency  amplitude.  (Greek  subscripts  refer  to  Fourier  components.  ) 
We  will  then  use  (1)  to  fi..d  the  amplitude  M^v  of  the  exp  ivt  component 
of  polarization.  Because  we  are  considering  only  stimulated  effects, 
v  is  a  linear  superposition  of  any  three  of  (±)thc  electric  field  fre¬ 
quencies.  We  write 


M. 

IV 


(2) 


where  a  is  the  usual  linear  polarizability  and  is  the  third  order 

nonlinear  polarizability  for  the  case  v  =  wp  +  +  w6  •  We  maY  now 

expand  (1)  in  terms  of  matrix  elements  m^b  of  the  dipole  moment  mj 
between  states  a  and  b(c,  d,  etc.  )  to  obtain 
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where  the  intermediate  states  b,  c,  and  d  are  summed  over  all  energy 
eigenstates,  and  the  initial  state  a  is  summed  with  its  statistical 
weight  factor  Pa.  Here  V2  =  u5  +  an<*  the  symbols  (ijkf,  abed), 
etc.  ,  mean  the  products  mia^mj^cm]tc^m^ ^a,  etc.  ,  of  matrix  ele¬ 
ments.  The  frequency  uab  is  the  energy  of  state  Ja>  minus  the  energy 
of  state  |b>,  etc.  Each  factor  in  the  denominators  may  be  thought  of 
as  having  a  term  -irab,  etc.,  representing  the  width  of  the  particular 
transition  in  that  factor. 

In  all  the  Raman  H2  amplifier  and  osc’llator  arrangements  with 
which  we  are  concerned,  two  of  the  three  electric  fields  in  the  non¬ 
linear  factor  of  (2)  are  very  nearly  components  of  a  linearly  polarized 
plane  wave  pump  field  whose  frequency  components  are  at  ±  vQ  with 
amplitudes  Eo  exp  (-  ikQ  •  R)  and  Eg  exp  (ikQ  •  R),  respectively;  kD 
is  the  (real)  pump  wave-vector  --nd  R  is  the  position  of  the  molecule. 
Various  processes,  including  coupled  or  pure  Stokes  and  anti-Stokes 
Raman  and  Rayleigh  scattering,  and  the  dc  and  ac  Kerr  effects,  are 
identifiable  as  arising  from  one  or  more  of  the  terms  in  (3)  with  certain 
intermediate  state  combinations.  We  concentrate  now  on  the  coupled 
Stokes -anti-Stokes  waves  with  which  we  must  deal  in  the  H2  Raman 
maser. 


D.  Coupled  Wave  Equations  for  Stimulated  Raman  Scattering 

We  consider  here  the  case  of  interest  to  the  H2  Raman  maser 
where  all  wave  polarizations  are  very  nearly  parallel,  and  therefore 
only  one  space  component  cixxxxPyS  of  the  nonlinear  polarizability 
tensor  is  involved.  We  will  study  only  the  case  where  the  pump  wave 
is  essentially  unattenuated  by  the  scattering  and  the  only  important 
nonlinear  terms  in  (2)  are  quadratic  in  the  pump  amplitude.  That  is, 
two  of  the  three  subscripts  f3y6  refer  to  one  of  the  two  conjugate  pump 
terms.  The  third  subscript  must  refer  either  to  a  first  Stokes  or  anti- 
Stokes  wave  in  the  Raman  effect. 

We  must  look  for  such  terms  in  (3)  which  produce  polarization 
at  the  first  Stokes  and  anti-Stokes  frequencies,  assuming  that  the 
initial  state  a  is  the  ground  state  of  the  Raman  transition  and  that  no 
dipole  matrix  elements  exist  except  for  electronic  transitions  (a3  is  the 
case  for  H2)-  The  most  important  of  such  terms  are  resonant  because 
the  middle  factor  in  a  denominator  nearly  vanishes  when  the  state  c 
is  the  upper  Raman  level  (1/2  *  ±wac)*  however,  we  will  find  that  non¬ 
resonant  contributions  to  Qxxxx  maY  also  be  important  to  stimulated 
Raman  characteristics. 

The  polarizability  terms  of  the  form  (3)  which  are  quadratic  in 
pump  amplitude  will  couple  a  Stokes  polarization  wave 
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Pg  exp  -  i(k  •  x  -  vt)  +  c.  c. 


(4a) 


to  an  anti -Stokes  wave 

P  exp  -  if  (2k  -  k*)  •  x  -  (2v  *  v*)t]  +  c.c. 

These  waves  are  accompanied  by  electric  fields 

exp  -  i(k  •  x  -  vt)  +  c.  c. 


and 


(4b) 


(4c) 


Ea  exp  -  i[(2kQ 


* 
k  ) 


*  -  (2 v 


*)t] 


+  c .  c . 


(4d) 


Equating  the  terms  in  Maxwell's  equations  that  have  the  same  space  and 
time  dependence  gives 


p  = 

(x  +  X°  j 

jE  +  X 

* 

E 

(5a) 
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^  ss  sj 

1  s  sa 

a 

II 

* 
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X  E  + 

fx  +  X° 

)E  * 

(5b) 

a 

as  s 

V  aa  a, 

/  a 

where  Xg  and  X-jf  are  the  normal  linear  susceptibilities  and  the  non¬ 
linear  susceptibilities  Xsa,  etc.,  are  equal  to  p<igaL8a.  Here  p  is 
the  number  density,  Lsa  are  local  field  correction  factor'  approxi¬ 
mately  related  to  the  Stokes  (anti-Stokes)  indices  action,  ns(na), 

by  LSa  =  (n8  +  1)^  (nl  +  1)^/81,  etc.,  and  the  polarizabilities  asa 
are  obtained  from  (3).  The  resonant  parts  of  these  polarizabilities  are 
seen  to  be  (for  the  'ase  of  parallel  polarizations  and  a  molecule  cer¬ 
tainly  in  the  lower  Raman  level  a) 
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(6d) 


where  D  is  the  resonant  denominator  u>ca  +  v  -  vQ  +irca-  It  has  often 
been  assumed  that  ass  =  aSa  =  aas  =  &aa-  This  is  seen  to  be  the  case 
in  the  limit  where  all  the  wave  frequencies  are  much  lower  than  the 
energy  differences  ^bC’  uda'  This  would  also  result  with  the 

Plazcek  model  of  coupling  via  a  polarizability  change  that  is  proportional 
to  the  amplitude  of  the  Raman  vibration.  However,  with  most  materials 
even  more  than  with  H£,  this  limit  is  not  reached.  In  order  to  estimate 
the  differences  of  these  coefficients  for  the  Qqj  (1)  Raman  transition 
of  H£  we  assume  that  the  intermediate  state  sums  in  (6)  are  imitated 
by  pretending  that  there  is  only  one  intermediate  level  b  whose  energy 
difference  from  the  initial  state,  u)^a,  is  90,  000  cm-1  (a  value  which, 
if  used  to  calculate  the  linear  polarizability,  gives  a  good  fit  for  optical 
dispersion).  The  Qqj  (0)  frequency  w,--,  is  4155.2  cm'*  and  the  ruby 
pump  frequency  vQ  is  14,402  cm""  .  Witn  these  values  we  find 
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a  :a  :a  :a 
ss  sa  aa  as 


1:1.  01:0.  83:0.  95  , 


(7) 


which  will  appreciably  alter  predictions  of  relative  anti-Stokes  intensi¬ 
ties  and  the  gain  near  the  index-matched  emission  angles  for  Hp#  at 
pressures  around  10  atm  or  more.  Our  own  numerical  analyses,  how¬ 
ever,  indicate  that  this  effect  does  not  predict  any  gain  anomalies  such 
as  those  that  we  and  others  have  observed. 

The  nonresonant  corrections  to  (6)  from  (3)  are  of  many  physical 
origins.  We  have  studied  ,the  various  types  of  nonresonant  terms  and 
found  that  the  "electronic"  contributions,  where  b,  c,  d  are  all  excited 
electronic  states,  probably  dominate  the  "nuclear"  terms,  in  which  c 
is  a  component  state  of  the  ground  electronic  level.  Terms  at  v  of  the 
latter  type  arise  from  the  tails  of  neighboring  Raman  transitions 
(especially  from  (0)),  and  Kerr  and  Rayleigh  type  effects. 

As  a  result  of  two  investigations,  we  believe  that  She  electronic 
contributions  to  the  Xsa  are  important  to  any  quantitative  understand¬ 
ing  of  stimulated  Raman  effects.  First,  knowing  the  locations  of  the 
important  lower  lying  singlet  electronic  levels  of  H2  from  published 
spectroscopic  data,  we  have  estimated  that  some  energy  denominators 
of  (3)  can  be  quite  small,  nearly  resonant,  for  the  wave  frequencies  of 
interest.  Hence,  the  sum  over  all  intermediate  electronic  states  need 
not  be  small  compared  with  the  resonant  Raman  terms  exhibited  in  (6). 
Second,  in  the  Landolt-Bornstein  tables  a  value  exists  for  the  dc  Kerr 
constant  of  liquid  H2  (19°K)  taken  from  some  unpublished  measure¬ 
ments,  equal  to  4  x  10*  9  esu.  We  have  worked  out  the  quantum  theory 
of  this  Kerr  constant,  assuming  that  it  arises  entirely  from  molecular 
reorientation,  and  found  a  predicted  value  several  orders  of  magnitude 
lower  than  this  reported  value.  This  means  that,  either  the  reported 
value  is  in  error,  or  that  it  arises  from  purely  electronic  nonlinearities. 
If  we  assume  the  latter  and  suppose  that  the  pump  field  E0  is  equally  as 
effective  in  producing  a  polarizability  change  at  the  Stokes  frequency 
6a88  as  the  dc  field  in  the  Kerr  measurement  (this  could  be  false  because 
of  changing  energy  denominators),  the  6aSs  is  of  the  order  of  the 
maximum  of  Reass  of  (6)  at  10^  atm  pressure.  We  will  need  to  assume 
for  our  -»urposes  that  this  occurs  rather  at  10  atm. 

The  effects  of  these  nonresonant  terms  in  the  X^j  are  difficult 
to  guess  in  the  absence  of  aay  knowledge  of  their  relative  signs  or 
sizes.  They  are  not  generally  equivalent  to  index  changes  at  the  two 
frequencies  because  the  cross  terms  Xsa,  Xas  produc  ed  gain  and  quite 
different  effects  from  index  changes.  However,  when  the  on-axis  Stokes 
waves  are  far  from  being  index-matched,  as  they  are  at  the  pressures  (20 
to  100  atm)  at  which  our  H2  Raman  lasers  are  commonly  operated,  it 
appears  likely  that  the  main  effect  will  be  to  produce  a  simple  index 
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change  6nv  which  is  proportional  to  the  pump  intensity.  Qualitative 
agreement  with  experiments  would  result  if  we  assumed  this  change  to 
be  at  a  pressure  Patm  (atmospheres)  and  pump  I0MW/cm2, 


&n  » 
v 


2x10 


-9 


P  *  I 
atm  o 


(8) 


In  order  to  obtain  an  idea  of  the  tendency  of  such  an  index  change 
to  self-focus,  note  that  under  the  typical  operating  conditions,  P  ^  = 
100  and  I  =  20,  this  change  in  index  would  cause  total  internal  reflec¬ 
tion  of  a  beam  at  3  mrad  to  a  plane  at  which  the  change  occurred  sud¬ 
denly.  The  effects  of  these  nonre sonant  terms  on  beam  coherence  may 
well  be  comparable  with  the  focusing  effects  arising  from  the  gain 
variation  across  the  pump  beam.  These  latter  effects  have  been  dis¬ 
cussed  by  Kogelnik  and  Yariv  and  were  considered  earlier  in  this 
report  in  connection  with  the  observation  of  mode  pinching  and  pinning. 


E.  Maxwell's  Equations  for  Coupled  Waves 

The  nearly  collinear  field  and  polarization,  Stokes-anti-Stokes 
waves  (4)  obeying  the  relations  (5)  are  solutions  of  Maxwell's  equations 
if 


[  -kV  +  v2(n2  +  4ir  X  )]  E  +  4ir  v2X  E  *  =  F 
1  s  ss,J  s  sa  a  s 


I  -k2c2  +  v2(n*  +  4irX  )]  E  *  +  4irv2X  E  =  F 

l  3L  dL  aa  '  ^  a  a  a  a  a 


aa 


a  as  s 


(9a) 

(9b) 


where  k-  =  2k  0  -  k  and  va  =  2v0  -  v  .  The  driving  forces  Fg  and 
Fa  include  externafsignals  and  noise  terms  which  are  correctly  given 
by  applying  the  quantum  Nyquist's  theorem  for  parametric  amplifiers. 
We  will  be  interested  mainly  in  the  gains  and  oscillator  thresholds,  so 
that  it  will  be  sufficient  to  study  the  homogeneous  form  of  (9)  (i.  e  ,  to 
solve  the  determinant  of  the  coefficients  of  the  LHS  of  (9)  for  possible 
values  of  k  and  v).  Results  of  this  are  summarized  below. 


F.  Steady  State 

The  steady  state  solutions  for  the  infinite  plane  wave  pump  are 
obtained  by  finding  the  complex  k  =  k'  +  ik"  which  satisfies  (9)  for  a 
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giveii  real  positive  v.  Those  solutions  for  k  which  have  the  largest 
imaginary  parts  k",  and  hence  the  most  gain,  will  be  most  effective 
in  amplifying  external  signals  and  will  have  the  lowest  thresholds  for 
oscillation  for  a  given  amount  of  feedback.  Fixing  the  directions  of 
k  '  and  k  "  as  well  as  v  and  putting  the  complex  determinant  for  (9) 
equal  to  zero  gives  the  magnitudes  k'  and  k"  of  k  1  and  k  Since  the 
results  are  insensitive  to  the  direction  of  k",  provided  that  k"  «  k' 
and  k  "  is  directed  near  the  forward  direction  as  we  have  supposed, 
we  sKall  take  k"  to  be  parallel  to  the  unit  vector  u  in  the  direction 
of  k\  That  i9,  we  will  look  for  solutions  of  the  form 


k  =  uk  =  u(k'  +  ik")  .  (10) 

^  ew 


If  we  define  the  dimensionless  deviation  of  k  from  its  normal 
value  by 


G  =  [kc/(n  v)]  -  1 

S 


(11) 


and  define  the  dimensionless  nonlinear  susceptibility  parameters 

Z.a(v>  s  2”xsa(v)/(".“a»  <12> 


and  so  forth  for  Zss>  etc.  ,  and  also  define  an  index  mismatch  param¬ 
eter  for  the  small  angle  9  between  kQ  and  k  by 


V  =  i  [v  ^n^  -  (2v  n  -  vn  )^  -  20^  n  v  n  v]  /(n  v  )^  ,  (13) 

2  1  a  a  o  o  s'  v  o  o  s  •*  a  a 


then  the  determinant  of  the  coefficients  of  (9)  gives  the  quadratic  equation 
for  G 


(G  -  Z  )  (AG  +  V  +  Z  )+Z  Z  =0  (14) 

ss'  aa  sa  as  ' 

where  we  have  neglected  terms  of  order  k"/k',  0'1 ,  and  G  (which  is 
of  order  Z)  smaller.  The  parameter  A  =  (2v0vn0ng  ”  nsv^^nava^ 
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is  nearly  v«g/v  foi  H^,  gas.  Recall  that  the  linear  indices  of  .  efraction 
at  tne  frequencies  v,  vQ,  and  v.  are  n8,  nQ,  and  n&.  The  so-called 
normal  Stokes  index-matched  emission  angle  Q  is  obtained  by  putting 
V  =  0.  The  normal  unmatched  Stokes  (or  anti -Stokes)  gain  is  seen  to 
be  G  =  Zss  (or  AG  ~  -v- Zaa). 

For  the  \  1)  Raman  transition 


V(0)  =  8.  89  x  10'7  P  -  0.429  02, 

atm 


(15) 


yielding  a  matched  Stokes  angle  Qq 


1.43  P 


1/2 


atm 


mi 


and 


A  =  0.  553  . 


(16) 


It  i.c  convenient  to  write 


Zss  *  [fK0c/<2n.w>l  i  r/(w  '  v0  +  v  ~  lF)  +  Z8S1  <17) 


where  p  is  the  "peak  gain  constant"  for  uncoupled  Stokes  waves  and  is 
generally  taken  to  be  1.  5  x  10“  3  cm/MW  for  the  pressure  broadenedt 


Q01  (l)line  of  half  width  2F  and  Raman  shift  w/(2irc)=  4155.2  cm 


-1 


The  nonre sonant  component  of 
(8)  ere  com  ~.t,  is 


ss 


is  represented  by  z  ,  which,  if 

J  & 
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ss 
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(18) 


are 


Although  at  present  the  relations  of  the  Z8a>  zsa,  etc.  ,  to  Zgg 
not  well  known,  an  idea  of  the  quality  of  the  actual  gain  solutions  can 
be  obtained  by  assuming  Z&„  =  £  =  Zga  =  Zag  =  Z  =  Z'  +  iZ",  a^ 

would  follow  from  the  Plazcek  model  plus  equal  nonresonant  terms. 
The  determinantal  equation  (14)  then  gives  for  the  solutions  G, 
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2AG.  -  -V  -  Z(1  -  A)  ±{  [V  +  Z(1  +  A)] ' 
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(19) 
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When  zgg  <  0.  1  the  maximum  gain  occurs  when  V  is  as  large 
as  possible,  that  is,  at  angles  as  f:.r  from  index  matching  as  possible. 
However,  when  z  >  i,the  maximum  gain  occurs  near  the  angle  0m 
where  V  =  -Z'(l  f  A),  if  that  is  possible,  i.e.  ,  for 


6  w  [2.  07  x  10'6  P  .  +  1. 04  Z'] 

1  atm 
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(20) 


For  the  case  of  a  large  nonresonant  susceptibility  (z  <2  1).  the 
largest  G"  evidently  occurs  near  a  frequence  where  Z'  i3  most  posi¬ 
tive  (i.  e.  ,  for  w  -  vQ  +  v  ~-r  where  v  is  on  the  low  frequency  side  of 
the  Stokes  line  peak).  The  nonresonant  term  zS8  is  positive  and, 
above  several  atmospheres  pressure,  would  contribute  the  major  portion 
of  Z’  if  (18)  is  a  good  approximation.  This  nonresonant  portion  Z^®- 
would  be  ~  1.  17  10”®  zggI0  2.  3  x  10”^  In  Patm  therefore 
0_  ~  Qq  for  IQ  <  10^  MW/cm'.  For  our  operating  pressures  and  if 
(18)  is  correct,  we  see  from  (19)  that  the  root  G  with  the  most  gain 
is  nearly 


G+  =  i  1.345Z'  -i  0.405Z"  ,  (21) 


which  is  greater  than  the  normal  peak  Stokes  gain  -iZgg  (vQ  -  w)  =  GQ 
if  zgg  >  1.  For  very  large  zs8  »  1,  (P atm  »  ^e  largest  G'j  is 
at  the  frequency  given  by  (  v  -  v  -  w)/r  .  34  and  is 

G  ’’/G  -  1.  345  z  +  0.372.  (22) 

+  o  ss 


We  have  seen  after  (20)  that  at  the  moderate  Ip  (<500  MW/cm  ) 
at  which  our  lasers  operate,  the  an^le  of  optimum  gain  Q  is  only 
slightly  below  the  index  matched  angle  $Q,  which  for  a  typical  operating 
pressure  of  50  atm  is  10  mrad.  Furthermore,  the  gain  falls  off  very 
rapidly  away  from  this  angle.  If  the  pump  beam  divergence  were  much 
less  than  one  might  suppose  that  the  R  am  an  wave  „o”ld  not  take 
advantage  ol  optimum  gain  for  more  than  a  certain  leng;h  before  the 
coupled  waves  would  diverge  from  each  other  in  space.  This  would  not 
be  true,  however,  if  the  nonlinear  index  (ReG)  of  the  coupled  Stokes - 
anti-Stokes  pair  were  large  enough  to  confine  the  wave  within  the  bounds 
of  the  pump  beam  by  dielec‘ric  waveguide  effect.  This  confinement 
would  occur  rather  suddenly  as  the  pump  power  increased  and  would 
bring  on  a  sudden  increase  in  the  single  pass  gain.  It  is  in  such  a  way 
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that  we  propose  that  the  sudden  onset  of  larg-  stimulated  Raman  outputs 
(as  pump  power  is  increased)  may  be  under*  ..d.  A  more  quantitative 
understanding  obviously  would  have  to  come  from  a  study  of  waves  of 
finite  cross-sectional  area  rather  than  of  the  coupled  infinite  plane 
waves  we  have  been  considering.  However,  even  without  the  complete 
theory,  one  can  cbtain  a  qualitative  picture.  The  increase  in  index 
6n  required  to  totally  reflect  a  wave  at  a  plane  by  a  small  angle  0m 
is  0m?2.  If  (8)  is  roughly  correct,  this  means  that,  for  2  x  10"9 
IQ  F>atm~'l-  Patm  x  (that  is,  at  Iq  ~  5  x  102  MW/cm2), a 

threshold  for  dielectric  guiding  of  the  high  gain  waves  should  appear. 

This  crude  estimate  is  near  enough  to  the  observed  thresholds  of 
—  40  M’V/cmZ  to  make  promising  a  further  investigation  of  the  hypothesis. 
One  might  say  that  since  in  the  oscillators  and  amplifiers  the  com¬ 
plications  of  this  extra  "parametric"  gain  and  dispersion  are  unwanted, 
perhaps  one  should  work  at  low  pressures,  ~  5  atm,  where  the  Raman 
gain  is  not  much  reduced  but  the  parametric  effects  would  be.  Unfor¬ 
tunately,  transient  effects  appear  to  become  important  at  pressures 
below  several  tens  of  atmospheres,  for  the  reasons  we  discuss  in  the 
next  section.  Therefore,  nonresonant  parametric  effects  may  well 
appear  in  a  practiced  system. 


G. 


Transient  Behavior 


To  analyze  the  transient  buildup  of  a  wave  packet  after  a  plant; 
wave  pump  of  constant  amplitude  has  been  turned  on,  we  write  the  initial 
wave  packet  as  a  superposition  of  plane  waves  of  the  form  (4)  with  dif¬ 
ferent  real  k.  Maxwell's  equations  (9)  are  then  satisfied  by  a  complex 
v(Jc)  =  v'(k)  +  iv"(k)  for  each  k.  The  wave  packet  at  a  time  t  later  is 
given  by  superposing  the  original  k  components  with  each  modified  by 
a  factor  exp  iv(k)t.  As  usual,  a  smooth  wave  packet  whose  spatial 
extent  is  longer  than  [the  linewidth  cm'*]'*  usually  travels  with  slowly 
changing  shape  with  a  group  velocity  whose  vector  components  are 
Sv'/Skj.  The  packet  grows  in  amplitude  as  exp  (-v")t,  a  negative 
v"(k)  indicating  gain. 


To  find  v  it  is  useful  to  express  it  in  terms  of  the  dimension¬ 
less  parameter  U: 


v  =  (kc/n  )  (1  +  U), 
s 


(23) 


lrom  which  Maxwell's  equations  (9)  for  the  couple  waves  give 

(U  +  zss)  (AU  -  v(0)  -  zaa)  +  zsa  Zas  =0  (24) 
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where  A  and  V  are  as  given  after  (14)  but  evaluated  at  v  =  ck/n8  , 
vQ  =  ck0/n£  ,  and  va  =  cka/n|  so  that  their  values  are  still  as  given 
in  (15)  and  (16)  for  gas  nea  Raman  shifts  and  at  not  too  high 
pressures  (so  that  n8  -  1  «  1).  From  (17)  we  see  that  we  may  now 
write 


Z  =  G  [  (UQ  +  B  -  i)"1  +  z  ] 
ss  ol  SSJ 


where 


B  =  (w  -  v  +  ck/n  >/r 

O  8 


(25) 


(26) 


is  a  dimensionless  measure  of  the  distance  off  resonance  of  the  k  being 
considered  and 


Q  h  ck/(n  T  ) 
8 


(27) 


is  one -half  what  would  usually  be  called  the  "Q"  of  the  Raman  transition. 
In  the  pressure  broadened  region  >  5), 


Q  «  107/P  . 

atm 


(28) 


for  the  Qqj  (1)  transition  at  300*K. 

It  is  evident  from  (25)  that  TJ  generally  satisfies  a  complex 
quartic  equation.  The  four  roots  correspond  to  the  normal  Stokes 
and  anti-Stokes  elect-cmagnetic  waves  and  two  "vibrational"  waves 
in  the  matter  when  coupling  is  turned  off.  In  the  special  case  of  the 
classical  Plazcek  picture  for  the  matter  vibration,  there  is  an  obvious 
cancellation  and  (24)  reduces  to  a  cubic  equation  in  U,  since  there  is 
only  one  possible  matter  wave  in  this  picture.  In  any  case  we  are 
interested  in  the  solution  U  which  has  the  most  negative  imaginary 
part  U". 

We  will  limit  our  considerations  here  to  the  case  where  V(0= 0) 
is  so  large  that  the  nearly  pure  Stokes  root  given  by  U  h  ZJS  =  0  has 
most  gain.  This  is  the  case  of  interest  when  the  index  matched  angle 
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(where  V  =  0)  is  large  enough  and  the  nonlinear  index  effects  are  small 
enough  that  we  have  an  ideal  Stokes  amplifier  or  oscillator  employing 
axial  waves  of  very  small  transverse  wave-vector  content.  We  will 
then  see  how  wave  packet  growth  rates  saturate,  an  effect  which  sets  a 
certain  minimum  Pa^m  a*  which  efficient  operation  is  possible.  To 
study  this  Stokes  root  we  will  not  write  the  zgg  term  explicitly,  but 
will  simply  modify  the  index  of  refraction  to  include  it  (ng  +  G0zgg-*-ng) 
because  in  this  case  it  has  only  the  effect  of  a  negligible  simple  index 
change.  The  Stokes  root  Ug  is  then,  from  (24)  and  (25), 

1/2 

2QU  =  (i  -  B)  {  [  1  +  4  QG  /(i  -  B)2]  -1}  (29) 

s  o 


which  at  small  gains  (4QGq  «  1)  gives 


Us  *  GQ/(i  -  B), 


(30) 


which  varies  with  ck/n_  just  as  does  the  normal  Raman  resonance  and 
peaks  at  B  =  0  where  the  time  gain  per  radian  is  GQ  just  as  was  the 
epatial  gain  per  radian.  We  may  break  up  the  solution  (29)  into  expres¬ 
sions  for  the  real  and  imaginary  parts  Ug'  and  Ug".  For  small  B 
(i.  e.  ,  near  line  center)  these  become 


-QU  "  =  l  (\Tl  +  40G  -  1)  (31) 

SO  u  o 


U  '  =  -BU  "/n/  1  +  4QG  .  (22) 

so  so  o  ' 


Thus  we  see  that  when  QG  >  1,  the  peak  gain  U  "  continues 
to  rise,  but  more  slowly,  as  ^1Gq/Q  rather  than  as  GQ.  A  plot  of  this 
result  We.s  given  in  Fig.  3,  Section  II,  showing  how  the  time  gain  constant 
saturates  with  increasing  pump  powers.  GqQ~  1  means  that  the  normal 
steady  state  gain  per  centimeter  is  of  the  order  of  the  Raman  linewidth 
(in  cm-*).  The  group  velocity  for  an  axial  wave  packet  v  is 


v  =  (c/n  )  (1  +  U1  +  QdU'/dB),  (33) 

g  s 
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and  from  (32)  this  becomes 


n  va/c  =  1  -  7  l  1  -  (1  +  4G  Q)'1/2], 

3  5  6  ° 


(34) 


showing  that  the  group  velocity  drops  to  half  its  normal  value  in  the 
limit  of  high  gain. 

The  apparent  saturation  of  gain  for  a  wavepacket  arises  when 
the  buildup  and  decay  of  field  at  a  molecule  is  fast  enough  that  its 
Fourier  components  lie  outside  the  amplifying  bandwidth  (2r  rps)  of 
i  —  molecule.  It  was  found  experimentally  that  Patm  had  to  be  greater 
than  20  atm  if  no  gain  savuration  was  to  be  evident. 
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V. 


SUMMARY 


Liquids  were  found  to  be  much  more  limited  than  Hz  gas  in  po¬ 
tential  performance  in  Raman  lasers  because  of  nonlinear  index  effects. 
Symmetric -molecule  liquids  are  the  best  in  this  regard.  In  these,  the 
origin  of  the  nonlinear  index  in  ill  cases  of  interest  was  found  to  be 
molecular  redistribution,  dominating  electronic  nonlinearities  by  at 
least  an  order  of  magnitude.  We  were  able  to  show  theoretically  (Ap¬ 
pendix  I)  that  simple  dc  Kerr  measurements  will  give  the  optical  non¬ 
linear  index  very  closely  if  one  assumes  that  the  difference  in  indices 
measured  in  the  Kerr  effect  is  3/2  the  optical  nonlinear  index  change. 
Available  Kerr  measurements  show  that  the  very  best  liquids  for  Raman 
lasers  cannot  handle  more  than  several  megawatts  without  the  formation 
of  beam  instabilities. 

The  pro  otype  H2  gas  Raman  oscillator -amplifier  we  have  con¬ 
structed  and  tested  demonstrated  brightness  gains  of  20  dB  and  absolute 
brightness  output  of  ~  10^  W-cm“2  sr"*  by  the  end  of  the  contract 
period.  The  potential  output  of  10 16  W-cm’2  sr*^  was  not  yet  obtained 
with  the  5  J  ruby  laser  pump  pulse  available,  probably  as  a  result  of 
noncritical  technical  difficulties.  A  possible  theoretical  model  which  we 
have  proposed  of  the  anomalous  behavior  of  the  H2  Raman  amplifier  ob¬ 
served  under  certain  operating  conditions  could  be  checked  with  standard 
nonlinear  index  measuring  techniques.  Wc  have  found  no  fundamental 
difficulties  that  would  prevent  extension  of  the  technique  to  achieve 
brightnesses  ~10^  W-cm-2  sr-i  or  higher.  We  feel  that  our  work 
shows  that  this  Raman  technique  would  produce  brightnesses  up  to  10^ 
for  fewer  dollars  per  W-cm‘2  sr-l  than  any  other  means  presently 
available.  Table  V  summarizes  the  brightness  values  presently  obtained 
and  those  which  could  possibly  be  obtained  with  our  existing  laser  systems. 
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TABLE  V 
Laser  Brightness 


System 

Present 
Performance, 
W-cm'^  sr"l  x  10^ 

Possible  Performance 
and  Conuitions, 
W-cm-2  sr_l  x  10^ 

Ruby  Giant  Pulse 
Oscillator 

0,  4 

2.0  Shorter  pulse 
and  more  com¬ 
plete  transverse 
mode  control 

Ruby  Giant  Pulse 
Oscillator  and 
Amplifier 

0.  025 

2.  5  Best  current 
technology  for 
crystals  and 
laser  design 

H2  Oscillator 

0.  08 

0.  4  Improved  output 
coupling  and 
pump  techniques 

H2  Oscillator 
plus 

H2  Amplifier 

0.  1 

30  TJear  diffraction 
limited,  strongly 
saturated  output 
from  our  device 

W 
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APPENDIX  I 

EFFECT  OF  MOLECULAR  REDISTRIBUTION  ON  THE 
NONLINEAR  REFRACTIVE  INDEX  OF  LIQUIDS 


(Reprinted  from  The  Physical  Reveiw,  Vol.  163,  205-206,  1967) 


Reprinted  from  THt  Physical  Review,  Vol.  163,  No.  1,  205-206.  5  November .  196/ 

Matal  U  W.  S.  A. 


Effect  of  Molecular  Redistributtat  on  flit  Rob* 
is««r  SifncUn  Index  ©I  liquids.  R.  W.  Hell- 
wakth  [Phys.  Re' .  152,  156  (1966)].  We  have 
discovered  an  important  error  in  the  sign  of  a  term 
calculated  in  the  Appendix.  As  a  result,  the  question 
raised  in  this  paper  as  to  die  accuracy  of  the  Kirk¬ 
wood  superposition  approximation  for  calculating 
certain  configurational  averages  for  liquids  can  now 
be  amwered:  The  approximation  gives  the  wrong 
sign  at  liquid  densities  for  the  nonlinear  free-energy 
term  and  also  for  other  quantities  whose  signs  are 
known  on  general  grounds.  Specifically,  in  Eq. 
(A17),  for  7i  the  factor  (P+2w*/3)  should  be 
(P-2**/3).  Also,  in  (A21)  the  factor  0.5066  should 
be  replaced  by  0.922  and  in  (A22)  the  factor  —0.117 
should  be  replaced  by  —0.467.  We  have  re-ex¬ 
amined  the  small  correction  K*  and  ascertained  tha  l 
its  sign  k  negative,  so  that  it  contributes,  though 
very  little,  to  die  worsening  of  the  failure  of  the 
.uperposition  approximation.  With  these  correc¬ 
tion,  die  superposition  approximation  gives,  in¬ 
stead  of  the  exkting  Eq.  (38), 

Un/iNE'W)  -  (f,+2e*/3)A— 2w*/3 

+w*(32+136Arl67A1)/(«#r) 
+1.13pr(l  +K«/4) .  (38) 

From  the  definition  (28b),  Utt  must  be  positive. 
For  symmetric  molecules  (A-C),  the  right-hand 
side  k  seen  to  fail  to  stay  positive  as  the  density  p 
iocreyies’ beyond  i.c/r,  where  r/8  k  the  molecular 
volume.  Norma!  liquid  densities  run  between  2/r 
and  3/r.  The  above  change  in  Un  requires  that  die 
expression  (41)  for  «t  be  altered  to 


«s- 


wpeffi 


*(1— 4wpa/3' 


* 


4i*pcf 

AH - (32+1366+ 167  A*) 

45r 

8 t»pW(A  5) 


+ 


f4.51o*/*J,  (41) 


which  should  also  always  be  positive.  The  first  two 


columns  of  Table  1,  which  were  calculated  with  the 
wrong  Eq.  (41),  have  no  significance,  except  to 
indicate  roughly  the  relative  magnitude  of  the 
factor  paf/inr1  (1  —  4»fcr/3)-4,  which  will  exist  even 
when  a  mere  accurate  method  than  the  super¬ 
position  approximation  k  found. 

We  have  looked  for  other  obvious  failures  of  the 
Kirkwood  approximation  and  have  found  several. 
Kirkwood  [J.  Chem.  Phy*.  4,  592  (1963)]  calcu¬ 
lated  corrections  to  the  Ciausius-Mossotti  formula 
which,  in  the  limit  of  small  polarizability  a,  would 
be  exact  for  all  densities  if  one  used  exact  correlation 
functions  instead  of  the  superposition  approxima¬ 
tion.  These  corrections  change  from  positive  to 
negative  as  the  density  increases  beyond  around 
2/r.  Thk  contradicts  the  demonstration  in  the 
present  paper  that  the  Ciausius-Mossotti  formula 
is  a  lower  bound  for  the  dielectric  constant  Another 
typical  failure  of  die  superposition  approximation 
occurs  in  the  calculation  of  the  positive  quantity 
<?.fjD..-'D,m'y)/N-  <2.P..*V^)‘  for  which  it 
gives  [with  Eq.  (37)]  64*V/(45pr)-2*VA  a 
negative  result  when  p>32/15r.  The  failure  cannot 
be  traced  to  the  fact  that  with  Eq.  (Al),  fp4tJV 
deviates  slightly  from  pi,  etc. 

The  following  typographical  errors  should  be 
noted.  A  minus  sign  should  be  inserted  before  one 
side  of  (3).  In  the  exponent  of  the  right-hand  tide 
of  Eq.  (24),- "pi— it"  should  read  “ut—u".  In  the 
second  line  of  the  first  paragraph,  second  column, 
of  p.  161,  "Jr"  should  read  “Jr”.  The  term  4((«,))» 
in  Eq.  (A19)  and  the  left-hand  side  of  Eq.  (A23) 
should  be  divided  by  ctE'*. 

The  conclusion  k  reached  that  there  k  no  con¬ 
sistent  usable  approximation  scheme  tor  computing 
the  configuration  averages  required  to  estimate  the 
contribution  of  molecular  redistribution  to  the  non¬ 
linear  index  *t  of  liquids.  However,  Eq.  (41),  which 
comes  from  the  Kirkwood  superposition  approxi¬ 
mation,  should  be  accurate  for  molecular  densities 
below  those  at  which  the  inconsistencies  discussed 
above  appear. 
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I.  INTRODUCTION 


Around  1S75,  Kerr  discovered  that  when  a  liquid  is 
immersed  in  a  uniform  electric  field,  its  optical  index  of  refrac¬ 
tion  for  light  polarized  parallel  to  this  field  changes  by  an  amount 
6n  if,  that  is  different  from  the  change  6  n^  for  light  polarized 
perpendicular  to  the  field.  *  The  effect  exists  either  if  the  strong 
applied  field  is  constant  in  time  (DC  Kerr  effect)  or  if  it  oscillates 
at  a  frequency  too  high  to  be  followed  by  any  nuclear  motions  (AC  Ken 
effect).  In  the  first  case  reorientation  of  permanent  dipole  moments 
may  contribute  to  the  effect,  but  i.ot  in  the  second  case.  We  shall 
refer  in  what  follows  only  to  cases  where  a  permanent  dipole  moment 
either  does  not  exist  or  plays  no  role.  Then  the  Kerr  effect  is 
generally  explained  as  arising  from  the  tendency  of  asymmetric 
molecules  to  lower  their  interaction  energy  with  the  strong  field 
by  reorienting  themselves  so  that  their  axes  of  maximum  polariza¬ 
bility  tend  to  be  along  the  field.  ^  However,  the  Kerr  effect  in 
liquids  of  nominally  symmetric  molecules  is  comparable  to  that  for 
molecules  possessing  an  asymmetric  polarizability  tensor.  This  is 
illustrated  in  Table  I  which  lists  the  size  of  the  DC  Kerr  effect  for 
some  typical  symmetric  and  asymmetric  molecules  in  terms  of  the 
parameter  Bo> 
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which  is  called  the  Kerr  constant.  Here  is  the  time  average  of 

the  square  of  the  strong  DC  or  AC  field  and  \  is  the  vacuum  vave  - 
leigth  of  the  "test"  beam  experiencing  an  index  change.  We  have 
pointed  out  previously  that  the  probable  origin  of  the  Kerr  effect  in 
symmetric -molecule  liquids  is  the  spatial  "redistribution"  or  "rear¬ 
rangement?1  of  the  molecules  under  the  influence  of  the  anisotropic 

3 

dipole-dipole  forces  induced  by  a  strong  applied  electric  field.  How¬ 
ever,  we  have  previously  attempted  to  calculate  only  the  magnitude  of 
the  quickly  responding  parallel  change  in  index  6  n ^  on  this  hypothe¬ 
sis,  ^  rather  than  the  much  more  easily  measured  difference  6  n jj  -  6n 

(i.  e.  ,  Kerr  effect). 

In  these  lectures  we  will  discuss  the  physical  basis  of 
molecular  redistribution  and  why  it  is  likely  to  dominate  other 
effects  which  might  also  contribute  to  the  Kerr  effect  in  all  but  a 
few  symmetric-n.  ecule  liquids  (Section  II).  We  then  develop 
general  expressions  for  ft  n  ^  and  ft  n^  in  terms  of  configuration 
averages  in  the  (zero  field)  equilibrium  ensemble  of  molecules 
(Section  III)  and  discover  that  the  result  depends  mainly  on  the 
equal-time  joint  probabilities  for  two,  three  and  four  particles 
to  be  in  various  close  configurations.  Th«s  Kerr  effect  is  probably 
the  most  sensitive  of  all  easily  measured  properties  of  a  fluid  of 
symmetric  molecules  to  this  statistical  aspect  of  short  range  order 
From  the  derived  expressions  we  show  that,  independently  of  the 
nature  of  these  correlation  functions,  the  ratio  6  n^  /  6  n  is  nearly 
-2  for  pure  molecular  redistribution  of  symmetric  molecules  (Sec¬ 
tion  IV).  It  was  exactly  -2  in  the  old  theories  for  asymmetric 
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molecules  tindergoing  no  redistribution.  Although  we  find  that 
the  ratio  is  not  -2  in  general,  we  show  that  in  any  normal  liquid 
the  ratio  will  probably  be  close  enough  to  -  2  so  that  one  may 
infer  the  value  of  the  fast  responding  part  of  the  nonlinear  index 
6n^(which  generally  governs  the  optical  self- focus  sing  properties 
of  the  liquid)  from  DC  Kerr  constant  measurements,  to  within  about 
20%.  This  is  important  since  other  ways  of  measuring  this 

accurately  are  much  more  difficult.  Finally  we  shall  outline  how 
Kirkwood's  "superposition  approximation"  may  be  used  to  evaluate 
the  required  averages.  We  find  that  the  result  of  this  calculation 
gives  the  wrong  sign  for  Bq  at  densities  about  roughly  half  liquid 
density  (Section  V).  The  results  are  given  for  molecules  of  as  sym¬ 
metric  polarizability  but  symmetric  (unperturbed)  intermole cular 

force  ),  and  are  seen  to  reduce  to  the  results  of  Debye  and  others 
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in  the  low  density  limit.  Unfortunately,  there  is,  at  present,  no 
way  known  to  calculate  the  required  configurational  averages  for 
liquids  in  a  self-consistent  manner. 
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II.  ORIGIN  OF  KERR  EFFECT  IN  SYMMETRIC- 
MOLECULE  FLUIDS 

On  the  basis  of  evidence  available  to  date,  it  would  appear 

that  symmetric -molecule  liquids,  such  as  those  of  Table  I,  whose 

-24  3 

polarizabilities  are  well  in  excess  of  10  cm  owe  their  Kerr 
effect  mainly  to  molecular  redistribution.  Therefore  the  remaining 
parts  of  these  lectures  will  be  devoted  to  a  study  of  this  mechanism. 
However,  for  liquids  such  as  of  helium  and  neon  whose  atoms  have 
low  polarizabilities  the  Kerr  effect  will  probably  be  mainly  due  to 
the  electronic  nonlinearity  of  the  atoms.  The  evidence  is  oriefly 
as  follows. 

Stuart  and  Yolkmann  observed  that  upon  increasing  the 

temperature  of  liquid  CCl^  14%  (to  339°K)  the  Kerr  constant  dropped 
5 

lOv'o.  This  behavior  is  inconsistent  with  electronic  nonlinearities 
being  dominant  and  is  consistent  with  what  is  known  from  redistribu¬ 
tion  theory.  On  the  other  hand,  Boyle,  et.  ah,  have  measured  the 
Kerr  constant  of  helium  gas  and  found  it  to  be  proportional  to  the 

density  as  would  be  expected  if  electronic  nonlinearities  were  domi  • 

6  7 

nant.  They  and  Langhoff,  et.  aL,  have  also  calculated  values  for 
the  electronic  nonlinearity  of  helium  and  predicted  values  close  to  each 
other  and  to  that  inferred  from  the  Kerr  measurements.  The  theory 
of  molecular  redistribution,  which  is  expected  to  be  ac  -urate  at  the 
densities  of  gas  employed  in  the  experiments,  predicts  a  nonlinearity 
that  would  be  much  smaller  than  that  observed  and  also  vary  as  the  square 
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of  the  density.  ^  Lunghoff,  et  al. ,  also  calculated  the  nonlinear 
polarisabilities  for  many  S-»tat*  atoms  and  ions  in  addition  to 
helium.  From  their  result#  and  a  crude  extrapolation  of  redistri¬ 
bution  effects  to  liquid  densities,  we  expect  (1)  that  the  Kerr  effects 
in  liquid  neon  would  also  be  electronic,  (?.)  that  in  liquid  argon  both 
electronic  and  redistribution  effects  m?  *ht  be  observed,  and  (3)  that 
in  more  polarizable  liquids  redistribution  is  likely  to  dominate.  For 

example,  if  we  use  the  calculated  nonlinear  polarizabilities  of  the 
4-4  -  >• 

C  and  Cl  ions  to  calculate  the  electronic  contribution  to  the 

Kerr  constant  JJL  of  C  Cl,  we  obtain  a  value  two  orders  of 
o  4 

magnitude  lower  than  what  is  observed  (and,  df  course,  do  not  pre¬ 
dict  the  observed  temperature  dependence  of  BQ). 

There  are  other  mechanisms  which  contribute  to  the  of 
liquids,  but  they  generally  contribute  a  negligible  amount.  Trie 
naturally  existing  different  isotopes  of  Cl  and  Br  would  cause  some 
XC14  or  YBr^  type  molecules  to  be  slightly  asymmetric.  This 
asymmetry  can  be  estimated  from  the  isotopic  shifts  in  spectral 
lines  and  appears  to  contribute  negligibly  to  B^.  This  is  not  surprising 
in  view  of  the  fact  that  a  molecule  like  cyclohexane  whose  polarizability 
and  number  density  in  liquid  are  very  close  to  those  c  f  C  Cl^  (4  4% 
nd  -10%  respectively),  but  whose  assymmetry  is  quite  large,  can 

nonetheless  exhibit  a  lower  Kerr  constant  than  C  Cl,. 

4 

Similarly,  the  electronic  asymmetry  of  those  molecules  that 
happen  to  be  thermally  excited  to  asymmetric  vibrational  states  can 
be  estimated  and  is  not  nearly  large  enough  to  account  for  <u  served 
Bq values  or  their  temperature  dependence. 
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Also,  if  the  light  whose  index  change  is  being  measured 
has  a  frequency  within  ~  100  Raman  (or  Brillouin)  ’’..e  widths  from 
the  frequency  of  the  strong  Kerr  field,  then  the  dispersion  from  the 
wing  of  the  Raman  line  may  affect  the  Kerr  constant.  This  "Raman 
wing"  will  never  significantly  affect  "DC"  Kerr  constants  but  could 
contribute  to  "AC"  Kerr  constants  and  might,  in  particular,  add 
significantly  to  the  nonlinear  index  seen  by  a  single  monochromatic 
wave.  Since  this  complication  is  easy  to  avoid  in  Kerr  measure¬ 
ments,  and  iis  easy  to  calculate  if  important,  we  will  not  consider  it 
further  here.  We  proceed  next  to  calculate  the  Kerr  constant  of  a 
fluid,  supposing  that  molecular  redistribution  were  the  only  mecha¬ 
nism. 
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HI.  RELATION  BETWEEN  A  KERR  CONSTANT  DUE 
TO  MOLECULAR  REDISTRIBUTION  AND 
SHORT  RANGE  ORDER  OF  A  FLUID 


Each  possible  way  of  arranging  N  molecules  inside  a 
volume  V  will  have  a  certain  probability  of  existing  (given  by 
statistical  mechanics)  and  will  exhibit  a  certain  wave  of  polariza¬ 
tion  density  Pj  exp(ijfe  -  iu>  t)  in  the  presence  of  a  propagating 

electromagnetic  wave  Ej  exp  (ijk  •  £  -  iut).  We  assume  that  each 
molecule  has  a  real  polarizability  a  that  is  strictly  linear  and 
isotropic  at  the  frequency  u  so  that  the  comp  'x  amplitude  m.1 
of  i.th  spatial  component  of  the  dipole  moment  of  the  y  th  molecule 
is  given  by  ^ 


(1) 


Y  Y 

The  local  electric  field  amplitude  E.  at  the  position  x  1  of  the 
yth  molecule  consists  of  the  transverse  propagating  electric  field 
plus  a  longitudinal  (Coulomb)  field  due  to  the  dipole  moments  of  all 
the  other  molecules: 


E 


.Us-*1 


Vij 


(2) 


where 
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(3) 


DijY  ^  -  6ij  I  3  *3(xiY  *  xiP)  'xjY  x/>  ■*,5|  5 

is  the  usual  dipole -dipole  matrix.  For  any  given  configuration 

V  Q 

{xT}  of  molecules,  \1)  and  (2)  yield 


m.V  =  uX..Y(3  E.e1^' 
i  ij  J 


where  w,ie  3N  x  3N  matrix  X^Y  is  defined  through  its  invei  se 


(X-1).  Y  P  =  6..6Y  P  +  aO.  Y 


The  polarization  density  amplitude  P.  for  the  (k,  w)  wave  is 
-1  Y  -ik*  xY 

V  c  ~  ~  and  we  wish  to  calculate  the  average  of  P.  over 

all  configurations.  If  we  express  this  average  as 


<  P>  =  x  •  -E- 

i  *ij  j 


then  the  susceptibility  tensor  X.j  is,  from  (4) 


xij  = 


where' 


Tij(  t  xY  }  )  =  X..Y  ^  (a/ V)  exp  [ ik-  (xY  -  xP)  ] 
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and  the  brackets  <>  indicate  either  a  quantum  or  classical  statis¬ 
tical  mechanical  average  over  configurations  {  x"^)  as  is  appropriate. 
The  average  Xjj  of  (7)  will  be  generally  a  function  of  the  existing 
electric  fields  because  the  intermolecular  induced-dipole- induced- 
dipole  interaction  energy  v  will  affect  the  weighting  function  used  in 
the  averaging.  Given  Xjj  it  is  a  trivial  matter  to  find  the  indices  of 
refraction  or  dielectric  tensor  of  the  liquid. 

•  1  1 

If  there  exists  a  strong  electric  field  E.  e*>*  *”w  '  in 

the  fluid,  whose  frequency  u  and  the  difference  ) -  w  |  are  large 

compared  to  any  intermolecular  collision  rates,  then  the  time  average 

o 

of  the  induced  electronic  dipole-dipole  interaction  energy  is 


1  i  i  *  i  v  . 

v  =  ■•  7  V  E.  E.  T..  (f  x*  }  ) 

4  l  j  ij  •  • 


(9) 


Y  ,  3  • 

for  any  configuration  {x  }.  The  prime  on  T.j  indicates  that  it  is 
the  same  operator  as  in  (8)  but  with  a  everywhere  replaced  by  the 

1  I  1 

real  polarizability  a  at  frequency  ta  and  k  replaced  by  k  . 

We  will  consider  here  only  classical  fluids.  Then  the  average 
<A([x^})>  of  some  operator  A  is 


v  +  v 


v+  v. 


(J*  dx1  •  •  •  dxN  Ae  }  /  (  Jdx1  •  •  •  d*Ne  kT  }. 


and  X..  may  be  expanded  in  v/kT  to  give  an  expression  exact  to 

J 

second  order  in  E*  : 
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(10) 


I  I  $ 

v..  =  <T..>  +  X--,  Ej  E 
*ij  lj  o  ljlm  x  m 


9.10 


ijlm  =  4kT  (<TijTlm  > 


-<  T.  >  <  T,  '>  ).  (11) 

o  ij  o  lm  o' 


The  average  <  >Q  indicates  the  average  over  the  unperturbed 
configurations,  where  only  the  unperturbed  intermolecular  potential 
v  acts. 

°  , 
For  small  k  and  k\  X  •  cannot  depend  on  k  or  k  . 

~  ~  ijlm  ~  ~ 

Therefore  the  rotational  symmetry  of  the  fluid  applied  to  (11)  with 
(8)  will  require 


X  -•  Y  -«  V  -4  V 

ljlm  ^jilm  *  ijml  jiml 


(12) 


or 


..S,  + 

lj  lm 


B(  6 ..  6  .  +6.  6  .,) 

'  ll  jm  lm  jl 


(13) 


as  k,  k  -  0,  and  the  two  coefficients  A  and  B  determine  x.-.  « 
w  w  ij  lm 

The  operator  T.j  is  too  difficult  to  evaluate  exactly  for  an 
arbitrary  distribution;  we  may  proceed  to  approximate  the  required 
averages  in  any  of  several  ways,  all  of  which  yield  the  same  result 
to  lowest  order  in  a.  We  could,  for  example,  use  the  fact  that  v 
may  be  expressed  as  a  minimum  of  a  certain  expression  and  mini¬ 
mize  an  approximate  expression  with  respect  to  some  parameters.  ^ 
Here  we  will  follow  the  program  used  by  Kirkwood  and  Yvon  to  calcu¬ 
late  corrections  to  the  Lor entz- Lorenz  formula  for  the  linear  diele- 
tric  constant  ot  a  liquid.  ^  The  RHS  of  (5)  is  rewritten  as  6  .  .6Y^(1  -  g  ) 

+  0  where  0^^  »  +  aD^,  ana  then  is  expanded 
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in  powers  of  0^  ^/(l  -  5). 


The  parameter  5  is  then  adjusted  to  be 


g  =  4irNo/(3V) 


(14) 


so  that  the  first  term  in  the  expansion  contributes  nothing  at  small 
k  to  the  transverse  part  of  the  linear  susceptibility.  The  hope  that 
this  produces  a  rapidly  converging  series  expansion  is  born  out  by 
the  high  accuracy  (~  1  to  5%)  of  the  Lor entz- Lorenz  formula  that 
then  results  from  the  0^^  =  0  approximation.  **  Expanding  (11) 
to  the  lowest  nonvanishing  order  gives  (p  =  N/V  is  the  average 
density) 


.  „  — 7j  «c..'^c'.  “,i>  -<c..^  <c',  "s ) 

ijlm  4kT(l  -{>*(!  1J  lm  0  *>  0  lm  0 


(15) 


where 


cijYP:  VPexpii‘ (16) 

and  similarly  for  C  but  with  k  replaced  by  k'.  The  next 

corrections  are  of  order  g,  g'  and  pT  smaller,  where  t  /8  ie 
the  molecular  volume.  Judging  from  the  accuracy  of  the  linear  term, 
we  feel  that,  if  evaluated  exactly,  (15)  would  give  the  nonlinear  part 
of  the  index  to  within  10%  for  typical  liquids. 
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IV.  6n  u  /  6nj^ 


The  most  simple  and  accurate  measurements  to  perform 
in  practice  are  of  the  Kerr  constant  Bq,  i.  e.  ,  of  in^  -in^  . 
However,  one  desires  the  value  of  6n^  alone  for  calculating  many- 
nonlinear  propagation  characteristics  such  as  the  self-focussing  of 
laser  light.  It  is  therefore  useful  to  note  that  to  within  the  accuracy 


of  the  approximation  (IS)  one  can  calculate  the  ratio  R=  6 n^ / 6 n^, 
and  hence  derive  fin^  from  BQ  (in  the  long  wavelength  limit),  with¬ 
out  knowing  how  to  perform  the  required  averages.  We  know  that  (lb)  must 
give  a  result  of  the  form  (13)  as  k-*0.  .But  EjC'j.  W>1=  ® 

so  that  multiplying  both  (IS)  and  (13)  by  and  summing  over  all 

space  indices  gives 


9A+  6B  w  0.  (17) 

I 

If  the  strong  electric  field  E.  is  applied  in  the  x- direction  (i  =  x) 
then 

R  .  Sn//  /  =  X^/x^  =  (A  +  ZB)/ A  (18) 

so  that  (17)  implies 

R  «  -2.  (19) 
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II  we  had  kept  nigher  ordez  terms  in  the  expansion  of  the 
exact  expression  (11).  A  would  not  be  exactly  -  2  B/3  nor  would 
R  be  -2.  But  one  might  expect  that  R  is  -2  to  within  10  or  20% 
for  even  the  densest,  most  polarizable  symmetric-molecule  liquids. 

We  will  see  from  the  discussion  of  the  next  Section  that  if  both  reorien¬ 
tation  and  redistribution  are  effective,  A  is  not  precisely  -2  B/3, 
even  in  lowest  order.  However,  again  it  is  close  enough  so  that, 
even  for  asymmetric  molecules,  R  would  never  be  likely  to  deviate 
from  -  2  more  than  10  or  20%.  It  is  well  known  that  when  only  reorien¬ 
tation  is  important  in  calculating  X^m  R  is  exactly  -2;  this  is  the 
famous  "Lange vin  relation".^  That  R  should  be  -  2  when  k  -•  0  implies 
that  the  trace  of  X^j  is  unchanged  by  the  application  of  the  strong  field. 
It  is  evident  that  this  should  be  true  when  there  is  only  the  reorientation 
of  molecules  whose  polarisability  tensor,  as  viewed  from  molecular 
axes,  remains  unchanged. 
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V.  EVALUATION  OF  X...  WITH  THE  KIRKWOOD 

ijim 

SUPERPOSITION  APPROXIMATION 

The  only  way  known  to  modern  statistical  mechanics  for 

estimating  the  averages  in  the  approximate  expression  (15)  for 

vijlm  **  to  U8e  t*le  Kirkwood  "superposition"  approximation  for 

the  two-,  three-  and  four-particle  correlation  functions  required 

there.  That  is,  one  must  take  the  probability  p^(1234)  of  finding 

12  3  4 

four  specific  particles  at  the  positions  £  ,  £  ,  £  and  £  (per 
fraction  of  the  total  volume  V  at  each  point)  as  a  product 

P^  (1234)  =  p  ( 1 2)  p(13)  p(34)  p(23)  p(24)  p<34)  /  Q  (20) 

over  all  possible  pairs  of  some  two  particle  function  p.  Then  the 
three-  and  two-particle  correlation  functions  are  determined  by  the 
identities 


p3(123)  =  Jp4(1234)dr4/V  (21) 


and 


P2(!2)  =J’p3(123)dr3/v 


(22) 


respectively.  The  normalization  factor  Q  ensures  that 
J,p2(12)dr1dr2  V-2  =  1. 


1  5 


The  only  known  function  p  for  which  dipole  averages  can 
be  carried  out  analytically  is  the  simple  step  function  of  the  absolute 
value  r  of  the  distance  between  the  two  space  points: 


p(r)  s=  1,  r  >  a 


=  0,  r  <  a. 


whence 


Q  =1  - 6t/N 


in  the  limit  of  large  N;  t  =4ira^  /3  is  eight  times  the  volume  of 
our  spherical  molecules. 

We  have  evaluated  the  averages  required  in  (15)  using  the 
above  correlation  functions.  We  have  also  carried  through  the 
averages  under  all  the  same  assumptions  except  that  the  molecular 
polarizabilities  are  tensors  (and  o  „  *)  that  depend  on  the 
orientation  of  the  yth  molecule.  This  means  that 
replaces  oD^^  in  the  definition  (5)  for  and  ^ik^°kj^ 

replaces  in  T„  and  the  interaction  energy  of  (9).  The 

average  polarizabilities  are  still  called  a  and  a'  and  the  "best" 

t  , 

5»  (  parameter*  are  still  4npa/3  and  4npa  /3. 

Also,  statistical  averages  must  now  include  the  appropriately 
normalized  integrals  over  all  orientations  of  all  molecules.  We  do 
not  however  change  the  unperturbed  distribution  functions  (20)  - 
(23)  to  reflect  any  molecular  asymmetry.  The  details  of 
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the  integration  techniques  and  the  organization  of  the  many  terms 
which  arise  are  given  in  reference  (3).  The  results  for  the  A  and 
B  coefficients  of  (13),  with  the  form  of  (15)  extended  as  above, 
are  (with  the  small  term  noted  in  ref.  3  omitted) 


A  - - - -j  {  A  (6  5'  5  -  y) 

4kT(l-?Hl-sr 


2  2 


1  r  (A  -  1)  8rr  p  on'  .  ,  •  3  128ir  poo'  „  13.  ,31.24l1 

2 1— 3 — B — “  +4.51oo  p  T  + -jg* (l"”4A+T6A  )J) 


(25a) 


and 


B  = 


- - Tjr  (M4  -  355') 

4kT(l-?Hi-5  ) 


4  51oa,p3T  t 


(25b) 


The  dimensionless  anisotropy  parameter  A  is,  as  usual, 


_2_  r  (a  -  b)  (a1  -  b1)  +  (b  -  c)  (b1  -  c1)  +  (c  -  a)  (c1  -  a1)  ■, 
~45l  ad1  J 


(26) 


where  a,  b,  c,  and  a',  b',  c'  are  the  polarizabilities  at  frequency 
(1)  and  w'  respectively  {.long  the  principle  axes  of  the  molecules' 
polarizability  tensor.  The  A  and  B  for  purely  symmetric  mole¬ 
cules  come  from  (25)  by  putting  A  =  0. 
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First  we  note  that  (25)  gives  the  nonlinear  index  coefficient 
we  have  calculated  previously  using  the  same  Kirkwood  approxima¬ 
tion  *  From  its  definition, 

n2  =  <  r>  v-xxxxfn  =  4  *( A+  2  B)/n  (2  ?) 

where  sill  primes  have  been  dropped  and  n  is  the  linear  index  of 
refraction  at  the  one  frequency  w  where  the  nonlinear  index  is  now 
self-induced. 

Secondly,  we  note  that  in  the  limit  of  low  densities,  only  the 
"A*1  terms  remain  and  we  have  the  classical  result  of  Debye  and 
others*  for  the  Kerr  constant  B  and  the  nonlinear  index  that 
arises  from  molecular  re- orientation  alone. 

However,  moV  important,  for  o  near  o',  X  becomes 
negative  for  pT  >  1.6.  For  liquid  densities  pt  runs  between  2  and 

2  2  i 

3.  Since  X  is  of  the  form  <T  >  -  <  T  >  when  o  =  o', 

xxxx  XX  O  XX  o 

we  see  that  x _  and  n,  ought  always  to  be  positive,  and  hence 

the  Kirkwood  approximation,  is  breaking  dot  a  entirely  at  liquid  densi¬ 
ties.  Other  failures  of  this  approximation  have  been  found  and  are 
noted  in  the  Erratum  of  reference  (3). 

The  formulae  (25)  give  a  maximum  positive  value  for 
:  mi  Bq  (at  p  T  0.  32)  which  is  typically  1/5  of  that  actually  obser  v<  i 
in  symmetric  molecules  at  liquid  densities.  We  see  therefore  that  at 
lost  densities,  before  (23)  fails,  the  calculation  is  givii.w  results  . 
oisteqt  with  what  one  might  expect  from  measurements  on  liquids. 
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Unfortunately,  measurements  of  BQ  of  symmetric  molecules  of 
high  polarizability  at  gaa  densities  where  (25)  should  be  accurate 
have  not  been  made. 

In  summary,  it  appears  that  molecular  redistribution  is 

the  main  cause  of  the  Kerr  constant  of  liquids  of  symmetric  mole- 

-24  3 

cules  whose  polarizabilities  are  much  larger  than  10  cm  . 
Experiments  at  different  densities  and  temperatures  wil-.  be  able 
to  isolate  the  contribution  of  molecular  redistribute  ?  to  the  Kerr 
constant  of  gases.  At  densities  lower  than  about  1/5  liquid  density 
(25)  should  give  this  contribution  accurately.  Since  this  contribu¬ 
tion  depends  only  on  the  short  range  order  in  the  fluid,  and  is  the 
most  sensitive  of  easily  measured  properties  to  the  short  range 
order,  future  theories  of  the  short  range  order  (particularly  of 
the  function  p^(1234))  may  be  readily  verified  by  Kerr  measure¬ 
ments. 
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TABLE  I 


Some  values  of  the  Kerr  constant  Bq  for  nonpolar  liquids 

in  units  of  B^*  the  Kerr  constant  of  benzene. 

_  o 

B^  =  4. 1  x  10  esu  for  \  =  546  mji . 


Liquid 

B  /B.  1 

O  D 

cs2 

8.6 

Benzene 

1 

Sn  Cl  4a 

0.  50b 

Si  Br4a 

0.43b 

cci4a 

0.  20b 

N2 

0. 19 

Cyclohexane 

0.  14 

Symmetric  molecules 

bThese  values  are  from  recent  measurements  of  N.  Georgf  and 
C.  Cook  (unpublished). 
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